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Introduction 

Recently, National Aeronautics and Space Administration (NASA) initiated a program to 
achieve the significant improvement in aviation safety. One of the technical challenges is the 
design and development of accelerated experiments that mimic critical damage cases 
encountered in engine components. The Nondestructive Evaluation (NDE) Group at the NASA 
Glenn Research Center (GRC) is currently addressing the goal concerning propulsion health 
management and the development of propulsion system specific technologies intended to detect 
potential failures prior to catastrophe. For this goal the unique disk spin simulation system was 
assembled at NASA GRC, which allows testing of rotors with the spinning speeds up to 10K 
RPM, and at the elevated temperature environment reaching 540°C (1000°F). It is anticipated 
that the facility can be employed for detection of Low Cycle Fatigue disk cracking and further 
High Cycle Fatigue blade vibration. The controlled crack growth studies at room and elevated 
temperatures can be conducted on the turbine wheels, and various NDE techniques can be 
integrated and assessed as in-situ damage monitoring tools. 

Critical rotating parts in advanced gas turbine engines such as turbine disks frequently operate at 
high temperature and stress for long periods of time. The integrity of these parts must be proven 
by non-destructive evaluation (NDE) during various machining steps ranging from forging bla nk 
to finished shape, and also during the systematic overhaul inspections. Conventional NDE 
methods, however, have unacceptable limits. Some of these techniques are time-consuming and 
inconvenient for service aircraft testing. Almost all of these techniques require that the vicinity 
of the damage is known in advance. These experimental techniques can provide only local 
information and no indication of the structural strength at a component and/or system level. The 
shortcomings of currently available NDE methods lead to the requirement of new damage 
detection techniques that can provide global information on the rotating components/system, 
and, in addition, they do not require direct human access to the operating system. 

During this period of research considerable effort was directed towards the further development 
of experimental facility and development of the vibration-based crack detection methodology for 
rotating disks and shafts. A collection of papers and reports were written to describe the results 
of this work. The attached captures that effort and represents the research output during the grant 
period. 
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ABSTRACT 

Rotor health monitoring and on-line damage detection are 
increasingly gaining the interest of manufacturers of aircraft 
engines. This is primarily due to the fact that there is a necessity 
for improved safety during operation as well as a need for lower 
maintenance costs. Applied techniques for damage detection and 
health monitoring of rotors are essential for engine safety, 
reliability and life prediction. 

This paper describes experimental techniques and 
preliminary analytical results related to crack detection in 
rotating disks. The concept behind the undertaken crack 
detection approach is based on the fact that the development of a 
crack results in a distorted strain field within the component. 
This in turn is associated with a minute unbalance. By 
conducting on-line monitoring of the disk’s radial vibration 
amplitude and phase, it is possible to see changes in the center 
of mass of the rotor system. 

To achieve the experimental set-up necessary to verify and 
study this crack detection technique, a unique disk spin 
simulation system was recently assembled by the nondestructive 
evaluation (NDE) group at NASA Glenn Research Center. The 
system allows for precision controlled spin tests that can 
facilitate the application of various sensing technologies for in- 
situ detection of rotor damage. A new, innovative capacitive 
sensing system was implemented to monitor blade tip clearance 
and the corresponding change in the center of mass of the rotor 
system. Two, separate 25.4 cm (10 in.) diameter disks were 
utilized for this study. The disks had machined teeth for 
imitating compressor or turbine blades. One was composed of 
titanium, while the second was nickel. In addition, the results of 
an FE analyses examined the behavior of healthy disks as well 
as cracked disks. 

The FE analyses were carried out utilizing various crack 
sizes and locations. The solutions focused on finding the 
changes in the maximum radial deflections as a function of 
rotational speed and crack characteristics. Utilizing the 
analytical results, the feasibility of measuring the changes of 
center of mass due to cracks in disks was defined for the given 
rotor system. 


INTRODUCTION 

Recently, the United States set the ambitious goal of 
reducing the fatal accident rate for commercial aviation by 80% 
within ten years [1]. In turn, the National Aeronautics and Space 
Administration (NASA) in collaboration with Federal Aviation 
Administration (FAA) and other federal agencies, universities, 
as well as airline and aircraft industries, responded by 
developing the Aviation Safety Program (AvSP). The AvSP 
provides research and technology products needed to help the 
aerospace industry achieve the challenge to improve aviation 
safety. The Nondestructive Evaluation Group at the NASA 
Glenn Research Center at Lewis Field in Cleveland, Ohio is 
currently addressing the development of propulsion system 
specific technologies intended to detect damage prior to 
catastrophe under the propulsion health management task. 

A unique disk spin simulation system was assembled at 
NASA Glenn Research Center. The system allows for disks up 
to 46 cm (18 in.) in diameter and rotational speeds in excess of 
10000 revolutions per minute (RPM). In addition, the system 
has elevated temperature capabilities of 540°C (1000°F). 
Various NDE techniques can be integrated and assessed as in- 
situ damage monitoring tools. Furthermore, controlled crack 
growth studies at room and elevated temperatures can be 
conducted on subscale spinning rotors. 

Currently, the focus of experimentation is the utilization of 
a non-contact capacitive method to measure the gap between the 
sensor and the rotating disk [2]. The sensor and the 
accompanying software are packaged as the HiBcind Vibration 
Monitoring System by Exsell Inc [2]. The sensor’s 
accompanying software calculates the blade tip clearance, the 
radial vibration, the blade time-of-arrival, the circumferential 
vibration, as well as the change in the center of mass of the rotor 
system. The change in the center of mass, calculated utilizing 
the displacement data, has been shown to be sensitive to cracks 
as small as 1.27 mm (0.05 in.) in jet engine rotors tested in spin 
pits [3], Although, at this point in time, the changes seen in this 
parameter have only been characterized in a subjective fashion. 
A deeper understanding of the relationship between crack 
progression and the change in the center of mass can be 
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achieved with the implementation of highly controlled crack 
initiation and growth tests on subscale spinning rotors. With the 
current disk spin simulation system such controlled tests can be 
conducted. In addition, other NDE sensing technologies like 
thermal acoustics and eddy current can be assessed. Before an 
experimental program can be initiated, an analytical modeling 
campaign must be undertaken. First, the current disks supplied 
with the spin system were analyzed. The disks were required to 
achieve rotational speeds of up to 10,000 RPM. In addition, the 
disks were analyzed concerning the in-plane compliance. There 
needed to be enough flexibility that if a crack exists in the disk, 
crack opening would occur as the rotational forces are increased. 
This crack opening would reveal itself as a shift in the center of 
mass of the rotor system or as a change in the maximum radial 
displacement. 

The goals of this study were first to assemble the disk spin 
simulation system. To address the concerns of rotor health 
monitoring, an innovative capacitive sensor was integrated to 
monitor the radial displacement of a spinning disk. Examples of 
the software and its success at monitoring damage in previous 
spin pit tests are shown in later sections. In addition, a finite 
element study was conducted to investigate the feasibility of the 
current disk design as well as the behavior when presented with 
a crack. In an effort to examine damage, various crack sizes at a 
given location were introduced into the disk model. The FE 
analyses were employed to track changes in the maximum radial 
deflections due to the various types of crack damage. Lastly, 
rotor dynamics analyses were conducted by integrating the shaft, 
bearings, and support structure. This was done to assure that the 
current experimental set-up simulates a real aviation turbine 
with operating speeds beyond the first critical. 

Lastly, it should be noted that the experimental set-up is 
rather unique in the arena of spin testing due to the fact that the 
horizontal shaft is supported by bearings on both ends of the 
rotor. This is a more realistic representation of an actual turbine 
as compared to a traditional spin pit. In spin pits, a vertical, 
hanging shaft is utilized with only the top having a bearing 
support. The rotor of interest is located between the top bearing 
and the bottom end of the shaft. Because the bottom end of the 
shaft is unconstrained, it is relatively free to move as a result of 
any dynamic forces. Therefore, it is assumed that any physical 
changes in the rotor will induce larger displacement variations in 
a spin pit set-up than in the more rigid, two bearing system 
employed here. 

EXPERIMENTAL METHODS AND PRELIMINARY SPIN 
PIT RESULTS 

The disk spin simulation system is shown in Figure 1. The 
major components of the system consist of an inverter driven 
induction motor (2.2 kW: 3 HP), double belt drive, rotating shaft 
with mandrel, ceramic element heater, and containment 
shielding. The motor pulley to shaft pulley ratio is 2:1. This 
allows for rotor speeds of up to 10,000 RPM. As mentioned 
above, the ceramic element heater can achieve rotor 
temperatures of 540°C (1000°F). 



Figure 1 Photograph displaying disk spin simulation system. Note that 
left side containment cover is removed for viewing rotor. 

A non-contact capacitive probe, shown in Figure 2, is 
utilized to measure the radial gap between the sensor and a 
rapidly moving grounded target. The rapidly moving object is 
the disk. The sensor is attached to the base of the spin system, 
approximately 2.03 mm (0.08 in.) from the tip of the blade. The 
innovative capacitive sensor is based on a DC offset rather than 
using modulation techniques. The DC voltage, in conjunction 
with the motion of the rotor, allows for bandwidths up to 50 
MHz [2]. These high acquisition rates are advantageous for 
monitoring individual blade vibrations in turbines, which in 
turn, assist in the quantification of high cycle fatigue in blades. 
By using an auto calibration technique, the sensor system 
eliminates most cable, vibration and temperature noise sources, 
thereby, allowing it to operate in harsh turbine environments. 
This auto calibration occurs during the nidi periods when no 
blades are in the vicinity of the probe. During this time, zero 
voltage is present. For a technical description of the capacitive 
sensor, the reader is advised to see reference [2]. 



j (electrically grounded) 


Figure 2 Photograph of capacitive sensor attached to the base of the 
spin system. 

The DC signal from the capacitive probe pre-amp is 
digitized using a 20 MHz, 12-bit A/D card. The voltage is 
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calibrated and then employed to calculate the blade tip 
clearance, the radial vibration, the blade time-of-arrival, the 
circumferential vibration, as well as the change in the center of 
mass of the rotor system. For this study, the focus was on the 
change in the center of mass calculation. As mentioned in the 
Introduction, this value shows promise as a rotor crack detection 
tool. The change in the center of mass is based on analyzing the 
magnitude and phase of the first harmonic, which is 
synchronous with the frequency of the disk rotation. The values 
for the magnitude and phase are obtained by applying a Fast 
Fourier Transform (FFT) to the blade displacement data. As 
mentioned above, when no blades are present the preamp 
produces zero voltage. Therefore, an absolute measure of the 
distance to the disk valleys between the teeth cannot be obtained 
because they are beyond the sensor’s range. The only real 
displacement data are at the blade tips where true voltages are 
induced in the sensor system. As a result, the FFT analysis is not 
applied to the raw data, but to the plots representing the 
individual gear teeth displacements. This will become more 
apparent in the discussion of the previous spin pit results on 
actual fan disks from turbine engines. 

Displayed in Figure 3 are unpublished results provided by 
Exsell, Inc. The figure provides a sampling of how the sensor 
and its software have been applied in the past. Shown are the 
data from a spin test of a large fan in an evacuated chamber 
using a multiple element capacitive sensor inside a spin pit. As a 
result of the multiple elements, each blade passing was 
represented by three pulses. The low cycle fatigue spin test was 
designed to induce cracks in an aircraft fan by cycling the part 
from top speed to a substantially lower speed, then continuously 
repeating the cycle. The sensor system was measuring and 
displaying the raw signal, blade tip clearance, time-of-arrival of 
the individual blades, and processing the vibration signal to 
produce plots that tracked the small changes in the balance of 
the system. Under the controlled conditions of the spin pit, the 
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Figure 3 Display for capacitive probe software. Example of data 
collected during a spin test of a fan rotor. The four plots on the right are 
data for a single rotation. The three peaks for each blade are due to the 
fact that the probe has three elements. The left side plots are based on 
the magnitude and phase of the first hannonic that is synchronous with 
the rotation of the rotor. 


change in the balance of the rotor system was assumed to be a 
function of low cycle fatigue cracks in the fan disk. 

The three graphs on the right hand side of Figure 3 monitor 
the change in the balance of the part. The change was later 
confirmed to be due to a crack that initiated in the disk. The 
unbalance of the disk was measured in real time. The change in 
the center of mass measurement (i.e., the change in the 
magnitude and phase of the first harmonic) displayed the 
difference between the initial balance of the rotating system and 
the current location of the center of mass. It is assumed that 
when cracks form, a distortion in the strain field occurs, in turn 
causing minute changes in the balance of the system. These 
changes are then displayed in the graphs on the right side of 
Figure 3. The top graph is a polar plot of the change in balance 
of the system, and the bottom two plots are the amplitude and 
phase of the difference in unbalance as compared to the baseline 
data. The polar plot simply shows the x and y components of the 
vector defined by the magnitude and phase of the first harmonic. 
A rise in amplitude and a stabilizing of the phase are typical 
indications of crack initiation and growth. This is assuming 
some internal movement of the rotating structure does not cause 
the unbalance. Again, for this particular case, the system was 
tracking a crack as it initiated and propagated within the disk. 
Because the data of Figure 3 has not yet been published, further 
details cannot be given at this time. 

At this point, the reader is reminded that the objective of 
this report was to assemble a disk spin simulation system. The 
description of the experimental methods provides an 
understanding for the directions taken during the analytical 
studies. The analytical results will allow for controlled 
experimentation concerning damage detection in cracked disks. 
Such results will be reported in future publications. 

FINITE ELEMENT ANALYSIS OF DISK WITH CRACK 

The finite element analyses were directed at understanding 
the behavior of the isolated disk as a function of crack size. Plots 
of the maximum radial displacement were constructed as 
functions of crack characteristics and rotational speed. Changes 
in the maximum displacement give an indication of the expected 
magnitude change in the center of mass. 

The 30-blade disks utilized in this study were made of CP 
titanium, grade 2, and nickel CM 400. For each disk, the outside 
diameter was 25.4 cm (10 in.); the bore and outside rim 
thickness was 0.953 cm (0.375 in.); the thickness of the web was 
0.318 cm (0.125 in.), and the cross-section and height 
dimensions of the blades were 0.318 cm x 0.953 cm (0.125 in. x 
0.375 in.) and 1.27 cm (0.50 in.), respectively. Lastly, eight 
threaded holes, 0.635 cm (0.25 in.) diameter each, were drilled 
through the disk; four holes were located at a radius of 10.2 cm 
(4 in.), and four additional holes were at a radius of 8.9 cm (3.5 
in.). The eight holes were spaced every 45°. The holes are 
utilized for possible mass attachments. The finite element 
analyses (FEA) of the rotating disk were conducted using the 
COSMOS FE software. In the analyses, parabolic tetrahedral 
elements were used with a non-dimensional global size of 
0.19462. 
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Figure 4 Geometry of the tested disk. 

For this study, a circumferentially oriented crack was 
selected and placed in the web portion of the disk. A crack 
oriented in this way is opened and propagated by radial stresses. 
This particular crack was chosen on the assumption that it would 
cause the most dramatic shift in radial displacement and, hence, 
change in center of mass. The crack tip was tapered and the non- 
dimensional crack element size was 0.097309. The radial 
distance from the disk center to the plane of the crack was 4.70 
cm (1.85 in.). This corresponds to the location of the maximum 
radial stress as obtained utilizing a plane stress solution [4]. 
Three different types of circumferential cracks were defined and 
investigated in order to find the resulting radial in-plane 
deformations at multiple rotational speeds (see Table 1). 


Table 1. Crack dimensions utilized in FE analysis. 


Crack Type 

Crack Geometry 

Circumferential 
Length (in.) 

Radial Width 
(in.) 

Depth 

Through 

Thickness 

(in.) 

I 

0.3 

1/64 

1/32 

II 

0.3 

1/64 

1/16 

III 

0.3 

1/64 

1/8 



The radial tip displacements of the healthy, un-cracked, 
disks are shown in Figure 5 for both the nickel and the titanium 
disks. These are the FE calculated maximum radial 
displacements due to the centrifugal expansion of the disk under 


rotational forces. As expected, the relatively stiff disks had 
rather small centrifugal expansions. 

Next, FE analyses were conducted utilizing the crack 
characteristics from Table 1. Figure 6 shows the difference in 
the maximum displacements between the cracked and un- 
cracked titanium disks. Since the behaviors of the nickel and 
titanium disks were very similar, only the titanium data is 
shown. The calculated values are extremely small, and cannot be 
experimentally measured with the available capacitive sensors. 
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Figure 6 Maximum relative displacement of un-cracked and cracked 
titanium disk. 

Based on the crack analyses concerning the given geometry 
of the disk, the selected size, the location, and the geometry of 
cracks, the types of cracks listed in Table 1 generate changes in 
the radial displacements that are very small and experimentally 
immeasurable. In order to simulate a more realistic scenario, a 
new disk imitating a turbine disk has been designed and is 
shown in Figure 7. The design was based on the stress results of 
the FE analysis, and optimized in terms of crack-induced 
potential in-plane deformations. Future work will involve a 
more in depth crack analysis as well as experimentation with 
cracked disks. 
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Figure 7 Optimized design of the simulated disk. 
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ROTORDYNAMIC ANALYSIS OF DISK SPIN 
SIMULATION SYSTEM 

The following rotordynamic analysis was conducted to 
further understand the dynamical behavior of the experimental 
set-up. As stated earlier, the system is meant to simulate a real 
turbine system. In most aviation based turbines, the system is 
expected to operate at some speed above the first critical. The 
disk spin simulation system consists of a shaft having a diameter 
of 1.27 cm (0.5 in.) and a length of 45.72 cm (18 in.), running 
on rolling element, deep-groove, single row bearings, which can 
be considered as isotropic with speed independent 
characteristics. The disk shown in Figure 4 (i.e., the original 
design) is mounted at the midspan of the shaft. 

The equations of motion for a rotor system are 

M'q + C(Q)q + Kq = F(t) (a) 

where q denotes the displacement vector, M is the symmetric 
mass matrix, the C(Q) is an asymmetric matrix which includes 
an antisymmetric, speed-dependent gyroscopic matrix, and K is 
the symmetric (in this case) stiffness matrix. The right hand side 
of Equation (a) represents the force excitation vector resulting 
from the system’s residual imbalance as well as any additional 
imbalance which might be induced by a crack. 

The system’s natural frequencies are the result of the 
solution of homogeneous version of Equation (a), and can be 
written as 

K = ~u i +jv i (b) 

or, in a single degree of freedom form as 

4 = — (c) 

vwf 

where a di = v t is the i' h damped natural frequency, and 

4) = , U ' ■ the corresponding modal damping ratio. The 

Ju 2 + v 2 

values of co di are directly used for construction of the Campbell 
diagram. The real part of A, t indicates the possible instability of 
the system. 



Figure 8 Campbell diagram for rotor test rig. 


For the rotor with the disk made of nickel alloy, the 
Campbell diagram is shown in Figure 8. The intersection points 
of synchronous excitation line IX with the natural frequencies 
loci determine critical speeds. 

The disk’s diametral and polar moments of inertia are 
29.209 and 58.34 lbm-in 2 , respectively. The mass of the nickel 
disk is 2.13 kg (4.7 lb m ). The first and second critical speeds are 
approximately 4070 RPM and 37575 RPM, respectively. For the 
system with the titanium-based disk, the corresponding critical 
speeds are 5461 RPM and 37787 RPM. 

• A I .F~QH .» I A 


Figure 9 Mode shape of vibration for the rotor with nickel disk at speed 
4070 RPM. 

The FE model of the rotor system with nickel disk with 
indicated stations and disk locations, and the overlapped 
vibration mode shape corresponding to the first critical speed, 
are illustrated in Fig. 9. These results show that for the current 
system, it is practical to operate in a region beyond the first 
critical. 

DISCUSSION AND CONCLUDING REMARKS 

A disk spin simulation system was assembled in order to 
deal with the goals of damage detection in aircraft engines. The 
two bearing system allows for subscale component testing at 
room and elevated temperatures. To address the concerns of 
rotor health monitoring, an innovative capacitive sensor was 
implemented as a tool for monitoring the radial displacement of 
a spinning disk. The radial displacement data can be utilized to 
calculate the change in the center of mass of the system. 
Preliminary data from spin pit tests of actual fan rotors have 
shown this to be feasible. 

Next, FE and rotordynamic analyses were conducted in 
order to fully understand the experimental set-up. Cracks of 
various sizes were introduced in the web portion of the current 
disk design. The FE results indicated that the selected 
circumferential cracks had little impact on the radial 
displacements of the disks at rotational speeds up to 15,000 
RPM. This would indicate that the changes in the center of mass 
due to cracks would probably be impossible to experimentally 
measure for the current simulated disk design and selected crack 
scenario. As a result a new disk was designed in order to 
optimize crack induced in-plane deformations. In addition, the 
new disk more closely resembles an actual turbine disk. 

The rotordynamic analysis was carried out to assure that the 
experimental set-up would operate in a speed regime beyond the 
first critical. The results demonstrated that the system is 
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properly designed and can traverse the first critical speed. 
Again, this simulates the operating regime of actual aircraft 
turbines. 

Future work in the area of rotor health monitoring will 
involve further analytical modeling of cracked rotors. As stated 
in the Introduction, the current methods of turbine rotor health 
monitoring are relatively subjective. Hence, there is a need to 
conduct an in depth modeling campaign to predict disk behavior 
in the presence of cracks. Such models will increase the 
understanding and reduce the subjectivity of this crack 
monitoring technique. The results of the models can then be 
compared to experimental crack growth data produced using 
subscale, spinning disks. In addition, the models can lead to the 
optimization of sensor location as well as assisting in the data 
analysis. 
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ABSTRACT 

Different approaches are used to sense and to localize a damage of rotating structures. Most of the methods take 
advantage of the dynamic behavior of the structural model [1-7]. This paper uses the modal and sensor norms, as defined 
in [8], to determine damage locations. The proposed approach allows localization of damaged elements within a 
structure, and provides information concerning the impact of the damage on the structure’s natural modes of vibration. 

Keywords: modal approach, damage detection, turbine disks. 


1. INTRODUCTION 

Current design practice for high-energy aircraft gas turbine rotors is based on an approach called “safe-life” method, 
which assumes that all material and manufacturing conditions that may influence the fatigue life of a rotor have been 
accounted for in prior laboratory specimen and full-scale testing, in addition, the final design is usually based on 
minimum properties. However, aircraft gas turbine industry experience has shown that the presence of certain material 
and manufacturing anomalies can potentially degrade the structural integrity of high-energy rotors. These anomalies 
occur infrequently and, therefore, are generally not present at critical locations during the testing of laboratory-based 
specimens or actual components. Undetected material and manufacturing anomalies represent a departure from the 
assumed nominal conditions and have resulted in several incidents, including the loss of a DC-10 at Sioux City in 1989 
[9]. 

Critical rotating parts in advanced gas turbine engines such as turbine disks frequently operate at high temperature and 
stress for long periods of time. The integrity of these parts must be proven by non-destructive evaluation (NDE) during 
various machining steps ranging from forging blank to finished shape, and also during the systematic overhaul 
inspections. Conventional NDE methods, however, have unacceptable limits. Some of these techniques are time- 
consuming and inconvenient for service aircraft testing. Almost all of these techniques require that the vicinity of the 
damage is known in advance. These experimental techniques can provide only local information and no indication of the 
structural strength at a component and/or system level. The shortcomings of currently available NDE methods lead to the 
requirement of new damage detection techniques that can provide global information on the rotating components/system, 
and, in addition, they do not require direct human access to the operating system. 

However, early damage (e.g., cracks) detection is not easily achieved during the operation of rotating machinery. The 
principal difficulty lies with the fact that incipient damage can propagate inside a part of machinery without producing 
detectable changes in its operational parameters. A crack initiating at a critical location within a rotating machinery 
component (e.g., shaft or disk) can be fatal for its operation. Typically, such a crack produces an undetectable change in 
the overall structural frequency response. For this reason, advanced failure detection methods are sought that can 
improve the detection of local structural changes in vital areas of rotating machinery [10]. 


1 Correspondence: Email: Jerzy.T.Sawicki(a), grc.nasa.gov : phone 1.216.433.2467; fax 1.216.977.7150 


The objective of this paper is to demonstrate the feasibility of detecting and locating flaws in rotating components 
utilizing the approach of modal norms and sensor norms [8], The proposed method allows for localizing the damaged 
portions of a rotating structure. It addition, it provides information about the impact of the damage on the natural modes 
of the damaged structure. 


2. MATHEMATICAL MODEL 

A rotating turbine disk is an example of disk-like structure which is cyclically symmetric in that its geometry can be 
defined by a sector which is repeated an integer number of times in the circumferential direction. For these periodic disk- 
like structures, no simple equation of motion or closed-form solution is available for the modal properties. However, this 
structure, as all others, can be approximated by a discrete, finite element, type of model. Such a structural model is 
characterized by its mass, stiffness, and damping matrices, as well as by the locations of the sensors and actuators. It is 
represented by the following second-order differential equation: 

Mq + Dq + Kq = Bu, (la) 

y = C q q + C v q (lb) 

In this equation q is the structural displacement vector of dimension n d ; u is the input vector of dimension r; y is the 
output vector of dimension s; and M, D, K are the mass, damping, and stiffness matrices, respectively, with dimensions 
n d x n d . The n d x r input matrix B, characterizes the actuator locations, the output displacement and rate matrices C 

and C v of dimensions sxn d characterize the displacement and rate sensor locations, respectively. The mass matrix is 
positive definite, and the stiffness and damping matrices are positive semidefinite; n d is the number of degrees of 
freedom, r is the number of actuators, and s is the number of sensors. 

The above equation can also be obtained in the modal coordinates using the modal transformation. For a small 
proportional damping let 0 J i be the ;th natural frequency and </> t be the ;th natural mode, or mode shape. Define the 

matrix of natural frequencies Q = diagico, ,co 2 ,...,oJ n ) , and the modal matrix ® = 0, ... <j> n ] , of dimensions n d x n , 

consists of n natural modes of a structure. A new variable, q m , called the modal displacement vector is introduced, such 
that q = <&q m . This transformation produces the modal mass matrix M m = ® r M<J> , and the modal stiffness 
matrix K m = <f> K <I> . Both matrices are diagonal. If the damping is proportional the modal damping matrix D m = <f> 7 /.)<t> 
is diagonal as well. 

Pre-multiplying Eq. (1) by <t> r , and subsequently by M~' gives the following modal model (see e.g., [1 1]) 

q m +2ZQ q m +Qrq m = B m u , (2a) 

y = C mq q m + C mv q m , (2b) 

The term Z = 0.5 M~'D m QT l is a diagonal matrix of the modal damping, B m is the modal input matrix, B m = B , 

while C mq = C tj <I> and C mv = C, <I> are the modal displacement and rate matrices, respectively. 

The modal equations (2a, b) can be re-written as a set of n independent equations for each modal displacement 

4,m + 2 CjCOfimj + cofq mi = b tni u , 

y, =c mq ,q m <+c mvt q mi . 


i = !,...,« 


(3a) 

(3b) 



where £7 is the /th diagonal entry ofZ. In the above equations y t is the system output due to the /th mode dynamics; 
b mi is the /th row of B m ; and c , and c mvi are the /th columns of C mq , and C mv , respectively. Next, c mi is defined as the 
equivalent output matrix of the /th mode 


- + c,„ 


CO, 


(4) 


The terms \\b mi and \\c mi are the input and output gains of the /th mode (see [7]), and p | denotes the Frobenius norm 

of x, i.e., ||x|| 2 = y]tr(x T x ) . In the case of multiple sensors c mi consists of ,v entries c .. , j=l,...,s, that correspond to each 
sensor. 


3. MODAL NORMS APPROACH 

System norms, such as H 2 , H,, or Hankel serve as a measure of system “size”. The basis of the undertaken approach is 
decomposition of the system norms in terms of its modal norms [8], In other words, norm of a MIMO (Multi Input Multi 
Output) system/structure can be represented as the root-mean-square (rms) sum of norms of a SISO (Single Input Single 
Output) system/structure. The H 2 norm of a system is used to characterize the system itself, along with its modes, and its 
sensors. Let G(co) be a transfer function of a system. The H 2 norm of the system is defined as 




(5) 


When utilizing modal coordinates, the equations of a flexible structure are uncoupled (see (3)). In turn, the norms of the 
/th mode and yth sensor can be approximated by the following equations: 




| M- lie 


mi ||2 || mi] || 2 


2VC 


CO, 


(6) 


Define the H 2 norm of the /th mode, G mi , as the root-mean-square sum over all sensors (see [7]): 


IIg ,|f=y g , 

II ™ In || v lb ’ 

j = 1 




Define also the H 2 norm of the /th sensor, G sj , as the root-mean-square sum over all modes: 


(7) 


l G ,|=ZK| j = c...,s. (8) 

i=l 

The H 2 norm of the total system is approximately the root-mean-square sum over either all its modal norms, or all its 
sensor norms: 
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The modal and sensor norms are used in the following damage detection procedure. 



4. DAMAGE DETECTORS 


Denote the yth sensor norm of a healthy structure by|G s/i/ , and the yth sensor norm of a damaged structure by|G J<? . 

The /tli sensor index of the structural damage is defined as a weighted difference between the /th sensor norm of a 
healthy and damaged structure, i.e., 



( 10 ) 


The sensor index reflects the impact of the structural damage on the /th sensor. 

Similarly, denote the /th mode norm of a healthy structure by||G mfo .|| 2 , and the /th mode norm of a damaged structure 
by | G mdj 1^ . The /th mode index of the structural damage is defined as a weighted difference between the /th mode norm 
of a healthy and damaged structure, i.e., 


G,„ 


O',,,, = J 
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( 11 ) 


The /th mode index reflects the impact of the structural damage on the /th mode. 

The question arises how to measure the sensor and modal norms. It follows from Eq.(6) that the norm is determined 
from the system’s natural frequencies, modal damping ratios, and the modal input and output gains. The gains, on the 
other hand, are determined from the natural mode shapes at the actuator and sensor locations. Thus the measurements of 
natural frequencies, modal damping, and the modal displacements at the locations of actuators and sensors for healthy 
and damaged structures allows for the localization of structural damage. Note that the natural mode shapes characterize 
the structure itself, since they are determined from the modal matrix <t> . Thus, they do not depend on the input time 
history, u (/), nor on the actual system deformations, q(t). 


5. EXAMPLE 


A disk of uniform thickness, which imitates a model of turbine disk, is rotating with constant angular velocity <x>. 
Stresses are generated in the disk due to centrifugal forces given by prca 2 , where p is the mass density per unit volume, 
and r is the radial distance from the center. Treating the centrifugal force as the body force in the radial direction, a finite 
element analysis can be conducted. As a numerical example, a disk with an inner diameter of 5.08xl0' 2 m (2 in.), an 
outer diameter of 0.508 m (20 in.), and uniform thickness of 0.635xl0' 2 m (0.25in.) is considered. The disk is rotating at 
(y=l,047 rad/sec (10,000 rpm). The material properties are as follows: £=210 GPa, G=79 GPa, v=0.28, and p = 7700 
Kg/m’. The finite element (FE) model consists of 64 shell elements, each having 6 degrees of freedom, and 80 nodes, as 
is shown in Fig. 1. The equivalent nodal force vector is calculated as ([12]): 


r h = 


f f N 

M 

J J 

W) 


Inrdrdz 


(12) 


where N is the shape function of the elements, and b r , b- are the body forces in radial and axial direction, respectively. 

The damage case is considered as a 20% reduction in element No. 9’s stiffness. Thirty four strain-gage sensors are 
placed within the first quarter of the disk. The body forces, described previously, are considered the excitation forces. 




Fig 1. Finite element model of the disk. 

The sensor and the modal indices are shown in Fig.2a,b. The sensor indices in Fig.2a indicate that the sensors s2, s5, si 1, 
and si 2, all located at the damaged element No 9, exhibited the largest changes. The modal indices in Fig. 2b show that 
the first mode was heavily affected by the damage. 




Mode Number 


Fig 2. Sensor and modal indices for the damaged disk. 


6. CONCLUSIONS AND FUTURE DIRECTIONS 


The illustrated modal norm approach is based on decomposition of the system norms in terms of its modal norms. The 
presented results show that the method is capable to locate the damage and to provide information about the impact of 
the damage on the natural modes of the damaged structure. The factor which might limit the application of this method is 
the requirement of a number of available sensors distributed over the structure. 

For practical utilization in real structures, extensive testing using valid assessment of the NDE techniques is necessary. 
Future work in the area of rotor health monitoring will involve further analytical modeling of damaged (cracked) rotors. 
Since the current methods of turbine rotor health monitoring are relatively subjective, there is a need to conduct an in 
depth modeling campaign to predict disks behavior in the presence of cracks. Such models will increase the 
understanding and reduce the subjectivity of crack monitoring techniques. The results of the models can then be 
compared to experimental crack growth data produced using subscale, spinning disks. In addition, the models can lead to 
the optimization of sensor location as well as assisting in the data analysis. 
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ABSTRACT 

This effort is being pursued in the Propulsion System Health Management 
(PSHM) area of the Single Aircraft Accident Prevention under the NASA Aviation 
Safety Program (AvSP). The goal of the AvSP is to reduce the aircraft accident rate 
by a factor of five within 10 years and a factor of 10 within 20 years. The PSHM 
goal is to eliminate propulsion system malfunctions as a primary or contributing 
factor to the cause of aircraft accidents. The purpose of this work is to establish 
/improve the tools for engine diagnostics and prognostics for measuring 
deformation and damage of rotating engine components at the ground level and for 
intermittent or continuous monitoring in engine environment. 

In this work nondestructive evaluation (NDE)- based technology is 
combined with model-dependent disk spin experimental simulation systems to 
monitor rotor damage in real time. In addition wireless eddy current is being 
developed for on- wing and just-in-time NDE engine inspection to provide deeper 
access and higher sensitivity that extends on-wing capabilities and improves 
inspection readiness. In the long run, these methods can establish the base for 
prognostic sensing while the engine is running without any overt action like 
inspection. This damage detection strategy includes experimentally acquired 
vibration- and capacitance-based displacement measurements, and analytically 
computed FEM- and modal norms-based models of well defined damages and 
critical mass imbalances in rotating disks/rotors. 
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BACKGROUND 


Uncontained engine failures and engine shut downs, relatively rare 
occurrences, can be very costly and may lead to tragic loss of lives. In the life 
management strategy and practices of commercial engine hardware, the safe life 
design philosophy dictates minimal inspections and “inspections of opportunity”, 
and requires specific inspections mainly in response to a field durability issue. 

Complete teardown of an engine can cost in excess of one million dollars. 
Hence, the capability to perfonn additional inspections on wing has both safety and 
economic benefits. Development of cost effective, on-wing inspection techniques 
for the detection of defonnation, fatigue damage, and cracking can have significant 
economic and safety benefits to the engine and airline industries and the general 
public through the prevention of propulsion system malfunctions. 

“Inspections of opportunity” includes florescent penetrant and advanced 
inspection techniques such as eddy current and ultrasonic scanning on engine parts 
when in the shop environment. In most cases, the defect of concern is fatigue 
cracking that may occur around complicated features and in difficult to access 
locations while in the on-wing configuration. To improve on- wing inspection 
effectiveness, requires development of both the sensor and the delivery mechanism 
to bring the probe to the inspection location and the proper modeling and simulation 
of the effect of defects for accurate remaining life estimates. 

Herein, newly established disk spin simulation systems and modal norms 
approach to modeling are described. Also preliminary results from wireless eddy 
current NDE technique are being reported on. 

HEALTH MONITORING OF ROTATING SUBSCALE ENGINE DISKS 

Two dual-bearing spin systems with horizontal shafts have been built at 
NASA Glenn for NDE-based and model-dependent disk spin simulation and 
analysis in order to improve the understanding of rotor/disk health monitoring by 1) 
mimicking damage in engines with actual subscale experiments, 2) optimizing the 
use of sensors type and locations, and 3) validating different modeling approaches. 




In figure 1 , one of the disk spin systems with high temperature capabilities 
is shown, and figure 2 shows typical results from capacitive type sensors, when 
monitoring tip blade displacements, due to center of mass changes initiated by 
deliberate mass imbalances that simulate crack initiation and propagation. 

Details on the disk spin systems can be found in [1]. 
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Figure 2. Display showing the results from mass imbalance tests 

More details on the modeling side can be found in [2] and in the following 
section. The NDE sensors employed in these facilities are capacitive DC type 
sensors for displacement measurements, eddy current and accelerometers for 
vibration measurements, and acoustic sensors for damage detection. 

MODAL NORMS APPROACH: AN EXAMPLE 

Different analytical approaches are used to sense and to localize damage in 
rotating structures. Most of these methods take advantage of the dynamic behavior 
of the structural model [3-9]. This work uses the modal and sensor norms, as 
defined in [10], to determine damage locations. The proposed approach allows 
localization of damaged elements within a structure, and provides information 
concerning the impact of the damage on the structure’s natural modes of vibration. 

A crack initiating at a critical location may produce an undetectable change 
in the overall structural frequency response. For this reason, advanced failure 
detection methods are sought that can improve the detection of local structural 
changes in vital areas of rotating machinery [1 1&12]. 

Since the current methods of turbine rotor health monitoring are relatively 
subjective, there is a need to conduct an in depth modeling campaign to predict 
disks behavior in the presence of cracks. Such models will increase the 
understanding and reduce the subjectivity of crack monitoring techniques. The 
results of the models can then be compared to experimental crack growth data 
produced using subscale spinning disks as discussed in the previous section. In 


addition, the models can lead to the optimization of sensor location as well as 
assisting in the data analysis. 

A disk of uniform thickness, which imitates a model of turbine disk, is 
rotating with constant angular velocity co. Stresses are generated in the disk due to 
centrifugal forces given by pro? , where p is the mass density per unit volume, and 
r is the radial distance from the center. Treating the centrifugal force as the body 
force in the radial direction, a finite element analysis can be conducted. As a 
numerical example, a disk with an inner diameter of 5.08x10' m (2 in.), an outer 
diameter of 0.508 m (20 in.), and uniform thickness of 0.635xl0' 2 m (0.25in.) is 
considered. The disk is rotating at 0=1,047 rad/sec (10,000 rpm). The material 
properties are as follows: £’=210 GPa, G=19 GPa, v=0.28, and /t= 7700 Kg/m 3 . The 
finite element (FE) model consists of 64 shell elements, each having 6 degrees of 
freedom, and 80 nodes, as is shown in Fig. 3. The equivalent nodal force vector is 
calculated as [13]: 


where N is the shape function of the elements, and b r , b z are the body forces in 
radial and axial direction, respectively. 

The damage case is considered as a 20% reduction in element No. 9’s 
stiffness. Thirty four strain-gage sensors are placed within the first quarter of the 
disk. The body forces are considered the excitation forces. 
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Fig 3. Finite element model of the disk 


The illustrated modal norm approach is based on decomposition of the 
system norms in terms of its modal norms. The presented results show that the 


method is capable to locate the damage and to provide information about the impact 
of the damage on the natural modes of the damaged structure (Figure 4). 




Fig 4. Sensor and modal indices for the damaged disk 


WIRELESS EDDY CURRENT FOR REMOTE CRACK DETECTION 

Coupling eddy current (EC) probe technology, with that of wireless 
communication, allows the insertion of the EC probe into existing access ports of 
jet engines. Using readily available borescope access ports, provides for internal 
inspections deep in the engine with un- tethered EC (Figure 7). 



Figure 7. Current eddy current practices showing access limitations and possible 
cable failures that need to be replaced with un-tethered access deep in the engine 


The power source for the EC probe can be a small "watch" size battery 
capable of powering the RF/microwave sub system for the time necessary for a 
measurement. A wireless probe module with the criteria of compactness, low noise, 
low cost, and compatibility with selected instrumentation is being fabricated. The 
Design considerations for the EC remote communication include immunity to 
jamming, power dissipation, frequency allocations, interface to detector probe, 
modulation format and reliability of transmission. The technology demonstration 
will include eddy current functionality as well as remote communication. The goal 
is for the eddy current inspection sensitivity to be comparable to a current similar 
‘wired’ inspection. 

Measurements were made in a PW2037 from port to port with and without 
engine rotation. Databases were taken from API to the engine front, API to AP2, 
AP2 to AP3, AP3 to AP4, AP9 to API 1 to the tail of the engine, where AP 
indicated “access port”. Transmission losses and “suckouf ’ regions where 
transmission is very low due to resonance in the structures were identified. Phase 
data were corrected. Typical magnitude plots as a function of frequency are shown 
in figures 8 and 9 for several AP configurations. 


Magnitude plot between ports in the 2-3 GHz range 



Figure 8. Magnitude plot between ports from 2 to 3 GHz 


Results indicate that the data can be acquired in frequency bands, which are 
reserved for industrial, scientific and medical purposes bands, called ISM bands. 
Also, it was found that since the frequency, is dependent on blade size, separation, 
and position, which varies throughout the engine, one frequency range is not 
sufficient for communication between all sets of ports. The data processor design is 


being studied and a prototype RF transmitting system for the remote probe is being 
built. 


Magnitude plot between ports in the 5-6 GHz range 



Figure 9. Magnitude plot between ports from 5 to 6 GHz 

Alert Service Bulletins (ASB), Airworthiness Directives (AD), and 
candidate parts with cracks have been identified for technology transition activities. 

CONCLUSION 

The approach to health monitoring of rotating subscale engine disks adapted at 
NASA Glenn was briefly described. In addition, modeling cracked disks with 
modal norms demonstrated that the method is capable to locate the damage and to 
qualify its impact from changes in the natural modes of the damaged structure. 
These rotor/disk health-monitoring facilities can now provide the opportunities to 
improve the understanding of field damage issues by experimentally simulating 
them with a whole campaign of modeling approaches like the center of mass 
changes, conventional rotor dynamics, and modal norms analysis. Further, wireless 
eddy current (EC) was found to be feasible for remote communication and on-wing 
inspection deep in the engine. It was also demonstrated, that communication data 
between ports can be acquired in frequency bands reserved for industrial, scientific, 
and medical purposes, called the ISM bands. More work is planned for wireless EC 
inspection sensitivity to be comparable to current “wired” inspection. 
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ABSTRACT 

The dynamic response of a cracked Jeffcott rotor passing through the critical speed with constant acceleration is 
investigated analytically and numerically. The nonlinear equations of motion are derived and include a simple hinge 
model for small cracks and Mayes’ modified function for deep cracks. The equations of motion are integrated in the 
rotating coordinate system. The angle between the crack centerline and the shaft vibration (whirl) vector is used to 
determine the closing and opening of the crack, allowing one to study the dynamic response with and without the rotor 
weight dominance. Vibration phase response is used as one of possible tools for detecting the existence of cracks. The 
results of parametric studies of the effect of crack depth, unbalance eccentricity orientation with respect to crack, and the 
rotor acceleration on the rotor’s response are presented. 

Keywords: crack, rotor, diagnosis, acceleration, resonance, unbalance, nonsynchronous response. 

1. INTRODUCTION 

In recent years, there has been an increasing interest for the development of rotating machinery shaft crack detection 
methodologies and on-line techniques 1 ’ 9 . Shaft crack problems present a significant safety and loss hazard in nearly 
every application of modern turbomachinery, particularly in the power generation industry, but not only. Recently, the 
Federal Aviation Administration (FAA), following several reports of shaft cracks, called 10 for additional NDE 
inspections of certain main rotor shafts in helicopters. Conventional NDE methods, however, have unacceptable limits. 
Some of these techniques are time-consuming and inconvenient for turbomachinery service testing. Almost all of these 
techniques require that the vicinity of the damage is known in advance, and they can provide only local information and 
no indication of the structural strength at a component and/or system level. Also, the effectiveness of these experimental 
techniques is affected by the high measurement noise levels existing in complex turbomachine structures. Therefore, the 
use of vibration monitoring along with vibration analysis have recently received increasing attention. 

An extensive review of the early literature on cracked shafts was published by Wauer 11 . Dimentberg 3 was apparently the 
first to report the effect of the rotating asymmetry on the shaft lateral vibration. Henry and Okah-Avae 4 presented the 
results of analog computer study and reported the subcritical resonance due to an interaction between the crack and the 
rotor’s weight. Mayes and Davies 5 studied the behavior of a cracked shaft model which took into account opening and 
closing of the crack as a stiffness step function. Later, they extended their work and suggested a method for calculation 
of change in stiffness due to a deep crack 6 . Gasch 741 modeled the breathing crack by a spring-loaded hinge, and 
numerically demonstrated a subharmonic resonance. Grabowski 9 used modal approach to the problem and demonstrated 
numerically strong dependence of the rotor vibrational behavior on the crack position. 

The previously published research results focused primarily on the study of cracked shafts subjected to constant angular 
speed. However, it has been well recognized 10 that vibration monitoring during startup or shutdown can reveal important 
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machine malfunctions, especially for machines such as aircraft engines, which are subjected to frequent starts and stops 
as well as high speeds and acceleration rates. The topic of transient cracked rotor response has been treated by only a 
few authors. Plaut et al. 13 analyzed the transient response of a simply supported, rotating, Eulei-Bernoulli shaft with a 
breathing transverse crack, during its passage through a critical speed at constant acceleration or deceleration. Recently 
Sekhar 14 " 15 investigated the transient vibration response of a cracked rotor passing through its critical speed, utilizing a 
simple hinge model for small cracks. He made an assumption that the vibrations remain small in comparison to the sag 
of the rotor under its own weight. If a cracked shaft rotates slowly under the load of its own weight, then the crack will 
open and close once per revolution. He noted oscillations which are developed near critical speed. 

In this paper, the Jeffcott rotor model is analyzed and a transverse crack is assumed at the middle of the shaft. The 
characteristic changes in the accelerating rotor phase response due to the crack and are presented as a possible diagnostic 
tool for crack detection in the accelerated rotors. The criterion for the opening and closing of the crack has been 
developed based on the angle between the crack centerline and the shaft vibration vector which makes possible one to 
ignore the rotor weight-dominance assumption. Stiffness weakening effects of cracked rotor in both strong and weak 
axes are taken into consideration for deep cracks . 


2. ROTOR AND CRACK MODEL 

The theoretical model, called the Jeffcott rotor , employs a flexible rotor composed of a centrally located unbalanced disk 
attached to a massless elastic shaft which is, in turn, mounted symmetrically on rigid bearings (see Fig. 1(a). The shaft 
does have a transverse crack running across its section and located close to the disk. The stiffness of the uncracked rotor 
system is symmetric (isotropic) and the damping due to the air resistance effect is assumed to be viscous. 



Figure 1 : (a) Model of Jeffcott rotor with crack, (b) shaft and crack position described in inertial and rotating coordinate systems. 

The origins of the inertial (x,v,z) and body-fixed rotating (^,r|,^) coordinate frames coincide with each other and are 
“attached” to the center line of the bearing supports. At any instant of time, the §-axis remains perpendicular to the face 
of the crack, causing that (^,T|,(;) frame rotates with the same velocity as the rotor. The crack position with respect to the 
vibration vector is determined by angle \|/, which changes continuously with the rotor’s speed. Finally, the rotor’s 
unbalance eccentricity £ is oriented at constant angle (3 with respect to the crack or i;-axis. For a rotor subjected to 


constant acceleration a, the spin angle changes, i.e.. 


<h(0 = <F(0) +oy +- 


a r 


where co (| is the initial speed (assumed to 


be zero) and <f>(0) = -(3 . 


The opening and closing or “breathing” of the crack, as the rotor rotates, results in a periodic change in the rotor 
stiffness. In order to not be restricted to the weight dominance assumption, where the dynamic displacement is smaller 
than the rotor’s sag, the angle \|/ (see Fig. 1(b)) is used to judge the opening and closing of the crack. This angle more 
practically describes the “breath” of the crack by taking into consideration the influence of the whirl speed on the closing 
and opening of crack. The rotor stiffness in the crack and cross directions can be written as 8 : 
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For linearly elastic crack and rotor deflections, the changes in the stiffness A and A may be assumed to be 
constant. For all but very large cracks, the stiffness change A K is small and often assumed to be zero. However, in this 


6 AK l 

study, for the case of the deep crack, its effect is included assuming that A K = — . The steering function /(t|/ ) is a 

6 

periodic function depending on both, the rotor spin and whirl speed, and is shown on Fig. 2 for small and deep cracks. 



Figure 2 : Crack model; (a) the closing behavior of (I) the hinge model for small cracks and (II) Mayes’ modified function for deep 
cracks, (b) opening and closing of the crack as a function of the shaft position along its orbit. 


For small cracks, the rectangular function for the hinge model 8 is used (see Fig. 2(a)), which describes the crack’s 
sudden opening and closing, by switching from 1 (open) to 0 (closed). This function can be approximated by the Fourier 
expansion as follows: 
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In the case of deep cracks, there is a smooth transition between the opening and closing of the crack and this is reflected 
by the Mayes’ modified function 6 , i.e., 


l + cos\|/ 
/<*) = — 5— 


(3) 


3. EQUATIONS OF MOTION 


The equations of motion for a Jeffcott rotor, with a cracked shaft, in a presence of gravity forces and unbalance 
excitation, and subjected to constant acceleration, can be expressed in inertial coordinate frame as follows: 
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where z and y are the displacements of the disk. The stiffness matrix K(q t ) (where q = (z yf ) is nonlinear, depending 
not only on time, but also on the position of the shaft center along the orbit. 

The above equations of motion in the rotating coordinate can be written as 
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where the transformation between the inertial and rotating coordinate frames takes the form: 
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Incorporating the hinge model for small cracks into Eq. (5) ( AK =0), normalizing displacements with respect to static 
deflection, and using nondimensional time yields: 
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where here f(\\i ) is a rectangular steering function (see Fig. 2(a)). 
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Similarly, one can incorporate into the equations of motion (Eq. (5)) the model for deep crack, assuming A = — --- 
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and Meyes modified steering function ;(\|/) = (see Fig. 2(a)), and write that 
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In the above equations the following definitions and nondimensional variables were employed (see Nomenclature 
section): 
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(see Fig. 1(b)). 


4. RESULTS AND DISCUSSION 

The presence of a crack on the shaft results in the stiffness variation as a function of time. In most rotor dynamic 
analysis, the system response is obtained in the form of steady state amplitude data at each operational speed. In the case 
of transient analysis, the dynamic response can be obtained using time marching techniques. In this study Runge-Kutta 
method was employed to numerically solve equations of motion (Eq. (7a) and (7b)). The nonlinear responses of rotor 

a 

passing through the critical speed with several values of constant angular acceleration ratios (y = — — ) and different 

CO" 



crack depths (A K ) were evaluated. For all cases it was assumed that damping ratio £ =0.05 and normalized unbalance 
e 

eccentricity is — = 0.05 . 

The change in the rotor phase response due to the acceleration and crack is presented in Fig. 4. In general, due to the 
crack the amplitude of the phase angle decreases for the given speed ratio and acceleration rate. In addition, it can be 
seen that the crack induces significant oscillations in the rotor phase esponse (see Fig. (4b)), especially for higher 
acceleration rates. Such a characteristic “saw-cut” pattern of the phase waveform can find potential application in the 
crack on-line detection monitoring systems. 




(a) no crack (b) A K = 0.4 

Figure 4: Phase angle of the accelerated rotor passing through the critical speed for y =0.01, 0.8x10 3 , and 10 4 , and (3 = 0. 


1 - 6 i [ 1 i 1 1 r 



04 L « • « ■ * 

0.2 0.4 0.6 0.8 1 1.2 

Speed Ratio (co/co n ) 

(a) no crack 



Figure 5 : Vibration amplitude in Z-direction of the accelerated rotor passing through the critical speed for y = 0.8x10' 3 and (3 = 0. 




Figure 6: Vibration amplitude of the accelerated rotor passing through the critical speed for y = 0.8xlCT 3 , (3 = 0 and A K = 0.32. 


Figures 5 and 6 show nondimensional vibration amplitude in Z-direction and the overall amplitude, respectively, for the 
crack depth of AK = 0.32. It is apparent that the shaft crack reduces the stiffness of the system, shifting the resonance to a 
lower speed. It can be noticed the presence of sub-critical (1/3, 1/2) response peaks and large increase of vibration 
response at the fundamental resonance peak when the crack exists. 
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Figure 7 : Response of accelerated rotor passing through the critical speed with crack and gravity for y = 0.8x10"’, (3 = 0 and AK = 
0.15, 0.28, and 0.37; (a) vibration amplitude in Z-direction, (b) vibration amplitude. 


The vibration amplitude in Z-direction and the overall vibration amplitude transient responses for different crack depths and 
constant acceleration rate are illustrated in Fig. 7. It can be noticed (see Fig. 7(b)) that Ihe whirl vector develops significant 
oscillations near the fundamental resonance speed. The zone of critical speed becomes wider and the vibrations level increases with 
the depth of the crack. As a result, for the deeply cracked rotors there is no exact critical speed. 



Presentation of the rotor transient response in rotating coordinates also provides useful information. Thus, Fig. 9 shows 
vibration amplitude in ^-direction of the accelerated rotor passing through the critical speed for crack depths AK =0.32 
and 0.36. In the absence of the crack the resonant vibration amplitude almost does not show-up. However, if there is a 
crack then the fundamental resonance is apparent and additionally the second and third harmonics are present, which is a 
crack-characteristic phenomenon. It can be seen that even small change in the crack depth causes that the rise of 
amplitude. 



Figure 9 : Vibration amplitude in ^-direction of the accelerated rotor passing through the critical speed for y = 0.8x10" 3 and (3 = 0; (a) 


no crack, (b) A K = 0.32, (c) AK = 0.36. 


The effect of unbalance eccentricity orientation on vibration amplitude of a cracked rotor is shown in Fig. 10. 
Orientation of unbalance affects the amplitude of the fundamental resonance, but does have a little influence on the sub- 
harmonics. The maximum resonant vibration amplitude exists when the eccentricity vector is along the crack centerline 

((3=0). When the unbalance is on the opposite side of the crack ( (3 =180°), then the effect of the crack is significantly 
diminished. In this situation, the resonance amplitude is only about 1/6 of the corresponding amplitude for the case when 

(3=0. 



Speed Ratio (to/co n ) 


Figure 10: The effect of eccentricity orientation on vibration amplitude of a cracked rotor for (3 =0°, 90"and 180°, 

AK = 0.35 , and y= 0.8xl0‘ 3 . 




5. CONCLUSIONS 


The presence of a crack in the shaft of the rotor system affects the dynamic response significantly. The changes of 

vibrational behavior can be used to predict and/or detect the crack in the shaft. The following conclusions can be drawn 

based on the results presented in this paper: 

1. The developed model enables to study the cracked rotor dynamic response with and without weight dominance. The 
model includes small and deep cracks, as well as cross stiffness effect (for deep cracks). 

2. Characteristics “saw-cut” pattern in the phase response waveform can find potential application in the on-line crack 
detection monitoring systems. 

3. The crack reduces stiffness of the system, shifting the resonance to lower speeds. This effect is magnified by the 
depth of the crack. It can he noticed the presence of subcritical (1/3, Vi) response peaks and large increase of 
fundamental vibration response when the crack exists. 

4. For the accelerating rotor the zone of critical speed becomes wider and the vibration level increases with the depth 
of the crack. The whirl vector develops oscillation near the critical speed. 

5. The orientation of unbalance eccentricity with respect to the crack centerline affects fundamental resonance and 
does have much less effect on subharmonics. 
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external damping coefficient 
gravitational acceleration 
uncracked shaft stiffness 
cracked shaft stiffness in ^ -direction 

cracked shaft stiffness in r| -direction 
mass of the disk 
position of bearing centers line 
shaft center 

inertial coordinate system; Z = z hsr Y =yf n„ 
constant angular acceleration 
angle between ^ and unbalance 

acceleration ratio; y = cx / oF 

stiffness change ratio ( = A K c IK) 

variation of stiffness in £, -direction caused by crack 

variation of stiffness in r| -direction caused by crack 

unbalance eccentricity of the disk 
external damping ratio 

static deflection of uncracked rotor, T| j( = g /or 

orientation of unbalance eccentricity with respect to andz-axis 

whirling angle corresponding to whirling speed 

rotating coordinate system; c, is in the crack direction 

angle between ^ and the line connecting bearing and shaft center 

angle between ^ and z axes 

angular velocity of rotation 

critical speed of uncracked rotor 
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Abstract 

This paper describes the analytical results concerning the detection of a crack in a rotating disk. The 
concept of the approach is based on the fact that the development of a disk crack results in a distorted 
strain field within the component. As a result, a minute deformation in the disk’ s geometry as well as a 
change in the system’s center of mass occurs. Finite element analyses were conducted concerning a 
notched disk in order to define the sensitivity of the method. The notch was used to simulate an actual 
crack and will be the method utilized for upcoming experiments. Various notch sizes were studied. The 
geometric deformations and shifts of center of mass were documented as a function of rotational speed. In 
addition, a rotordynamic analysis of a 2-bearing, disk and shaft system was conducted. The overall 
response of the system was required in order to design the experimental system for operation beyond the 
first critical. The results of the FE analyses of the disk indicated that the overall changes in the disk’s 
geometry and center of mass were rather small. The difference between the maximum centrifugal radial 
displacements between the undamaged and damaged disks at 8000 RPM was 0.00014 in. for a 0.963 in. 
notch length. The shift in center of mass was also of this magnitude. The next step involves running 
experiments to verify the analysis. 


Introduction 

Rotor health monitoring and on-line damage detection are increasingly gaining the interest of 
manufacturers of aircraft engines. This is primarily due to the fact that there is a necessity for improved 
safety during operation as well as a need for lower maintenance costs. Applied techniques for damage 
detection and health monitoring of rotors are essential for engine safety, reliability and life prediction. 
Recently, the United States set the ambitious goal of reducing the fatal accident rate for commercial 
aviation by 80% within ten years [1], In turn, the National Aeronautics and Space Administration 
(NASA) in collaboration with Federal Aviation Administration (FAA) and other federal agencies, 
universities, as well as airline and aircraft industries, responded by developing the Aviation Safety 
Program (AvSP). The AvSP provides research and technology products needed to help the aerospace 
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industry achieve the challenge to improve aviation safety. The Nondestructive Evaluation (NDE) Group 
at the NASA Glenn Research Center at Lewis Field in Cleveland, Ohio is currently addressing the 
development of propulsion system specific technologies intended to detect damage prior to catastrophe 
under the propulsion health management task. 

Currently, the NDE group is assessing the feasibility of utilizing real-time vibration data for detecting 
cracks in turbine disks. The data is obtained from radial blade tip clearance and shaft clearance 
measurements using capacitive or eddy current probes. The concept of the approach is based on the fact 
that the development of a disk crack results in a distorted strain field within the component. This, in turn, 
causes a small deformation in the disk’s geometry as well as a possible change in the system’s center of 
mass. The geometric change and the center of mass shift are indirectly characterized by monitoring the 
amplitude and phase of the first harmonic (i.e., the lx component) of the vibration data. Spin pit 
experiments and full scale engine tests have been conducted while monitoring for crack growth using the 
above crack detection methodology [2, 3]. Even so, published data is extremely limited and the basic 
foundation of the methodology has not been fully studied. This foundation should include a theoretical 
basis and experimental data to support the models. In addition, the few tests that claimed success have 
been limited mostly to the spin pit facilities. In spin pits, a vertical, hanging shaft is utilized with only the 
top having a bearing support. The disk or disks of interest are usually located near the bottom end of the 
shaft. Because the bottom end of the shaft is unconstrained, it is relatively free to move as a result of any 
dynamic forces. Therefore, it is assumed that any physical changes in the rotor will induce larger 
displacement variations in a spin pit set-up than in a more constrained, two bearing system. The analytical 
model studied here, and the upcoming related experiments, deal with a horizontal shaft supported with 
end bearings and a centrally located disk. It should be noted that this system is somewhat unique in the 
arena of spin testing due to the fact that the horizontal shaft is supported by bearings on both ends of the 
rotor. This is a more realistic representation of an actual turbine as compared to a traditional spin pit. 

As mentioned above, the crack detection methodology has only been characterized in a subjective 
fashion. In addition, the methodology is lacking an in depth modeling campaign. A deeper understanding 
of the relationship between damage progression and the change in the rotor’s center of mass due to the 
geometric modifications of the damaged disk can be achieved with the implementation of subscale rotor 
tests with controlled damage levels. The objectives of this study included the design of an optimal 
subscale disk to represent a full scale turbine disk; finite element analyses of undamaged and damaged 
disks in order to define the disk’s deformation and the resulting shift in center of mass; and rotordynamic 
modeling of the complete disk and shaft assembly to confirm operation beyond the first critical. The 
analytical results will be verified in an upcoming experimental study focusing on the design of the 
laboratory and the implementation of unique and innovative sensors along with the accompanying 
acquisition/analysis software packages. 


Finite Element Analysis: Disk Design and Damage Behavior 

The finite element analyses had two goals. The first goal was to design a disk that will achieve the 
experimental objectives of flaw detection and at the same time be representative of an actual bladed, 
turbine disk. The second goal was to understand the behavior of the isolated disk as a function of damage 
state. For this study, the damage was in the form of a notch to be placed in a high stress area of the disk. 
Such a location allows for maximum geometric deformation. Plots of the maximum radial displacement 
and shift in center of mass while imposing centrifugal loads were constructed as functions of notch 
characteristics and rotational speed. Notch tip stresses were documented and will be utilized to define the 
load limits during the experimental stage of the continuing study. 
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Figure 1 . — Geometry of subscale disk. 


TABLE I. — PROPERTIES FOR HAYNES NICKEL ALLOY X-750 


Density, lb/in 3 

0.3 

Modulus of Elasticity, ksi 

31000 

Poison’s ratio 

0.31 

Shear Modulus, ksi 

11830 

Ultimate tensile strength, ksi 

145 

Yield strength, ksi 

95 


Figure 1 shows the final design of the subscale disk. The design was based on the fact that the disk 
should safely handle rotational speeds up to 25000 revolutions per minute (RPM) in the undamaged state 
as defined by the materials yield stress. A thinned web area and gear teeth were implemented to imitate 
the web area and blades of a turbine disk. Note that the gear teeth were also needed due to the fact that the 
capacitive displacement sensors to be utilized during experimentation were designed to monitor radial 
blade tip clearance. The nickel alloy, Haynes X-750, was the material chosen for the disk. The material 
properties for Haynes X-750 are listed in Table I. The disk’s final design had an outside diameter of 

9.25 in. and a bore thickness and outside rim thickness of 1.00 in. and 1.25 in., respectively. The thinnest 
portion of the web was 0.10 in., with the cross-section and height dimensions of the blades being 

1.25 in. X 0.13 in. and 0.33 in., respectively. Lastly, eight holes, 0.20 in. diameter each, were placed 
midway in the rim section. The eight holes were spaced every 45°. The holes are to be utilized for 
possible mass attachments or used for future notch initiation points. During the design and damage study, 
the finite element analyses (FEA) of the disk were conducted using commercially available FE software. 
In the analyses, three dimensional brick elements were utilized with an auto-mesh feature to create the 
model. 

Upon completing the subscale disk design, a parametric analysis of the damaged disk was conducted. 
The damage was in the form of a circumferential notch located in the web region as shown in figure 2. 
The notch was utilized to represent a crack and is the damage parameter to be implemented in the future 
experiments. Note that a notch oriented in this way is opened by radial stresses. The circumferential notch 
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was chosen on the assumption that it would cause the most dramatic change in radial displacement and 
change in center of mass. The notch width of 0.15 in. was based on the wire thickness and burn area of 
the electrical discharge machining (EDM) wire. 

Figure 2 and Table II list the dimensions associated with the various notch lengths studied. For all 
cases, the radial distance from the disk center to the plane of the notch was 2.12 in. This corresponded to 
the location of the maximum radial stress as obtained utilizing the FE stress analysis of the undamaged 
disk. Figure 3 shows the stress state at the notch tip as a function of notch length and the disk’s rotational 
speed in terms of RPM. The results of the figure will be utilized in defining the experimental notch length 
and disk speed required to maintain a safe test environment. 


TABLE II. — NOTCH DIMENSIONS UTILIZED FOR FINITE ELEMENT ANALYSIS 


Angle, degrees 

6 

12 

18 

24 

26 

28 

32 

Length, in. 

0.222 

0.444 

0.666 

0.888 

0.962 

1.036 

1.184 



Notch 


Notch location: 
R = 2.1 2 in. 



Figure 2. — Notch dimensions for particular case of 1.184 in. notch length (see Table I). 



Figure 3. — Von Mises stresses at notch tip as a function of notch 
length and rpm. 
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Next, the focus was on the global influence of the notch on the disk’s overall displacement behavior 
while subjected to centrifugal loads. In particular, results were documented concerning the change in 
maximum centrifugal expansion as measured at the tips of the gear teeth, and the relative shift in center of 
mass when in the deformed state. These values correspond to values that will be experimentally measured 
using capacitive displacement probes that monitor the radial blade tip clearance and eddy current probes 
that monitor radial shaft displacements adjacent to the disk. Figure 4 displays the maximum radial 
displacement due to centrifugal expansion as a function of notch length and rotational speed, while 
figure 5 displays the net change after subtracting the displacement values of an undamaged disk. 


1.2x10-3 krpm 



Figure 4. — Maximum radial displacements as a function of notch 
length and disk rpm. 


3.5x10-4 



Figure 5. — Net radial displacements after subtracting displacement 
values from undamaged disk. 
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3.0x10-4 



-0.5 ' 1 1 1 * 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Notch circumferential length, in. 

Figure 6. — Radial shift of center of mass as a function of notch 
length and disk rpm. 


Viewing figures 4 and 5, it was apparent that the notch induced changes in maximum displacement 
were rather small. For example, the net change in the maximum radial displacement of a notched disk 
(notch length = 0.962 in.) at 8000 RPM was 0.00014 in. This particular set of parameters was of interest 
due to the fact that it will be experimentally feasible to test these conditions since the notch tip stress was 
maintained below the yield strength of the material (see figure 3). 

Next, parameters from the FE model (e.g., the mass matrix and node displacements while in the 
deformed state) were utilized together with the equations of basic engineering mechanics to calculate the 
disk’s shift in center of mass [4]. Figure 6 shows the approximate shift in the center of mass as a function 
of notch length and rotational speed. The values were similar in scale to the net change in radial 
displacement (see figure 5). Therefore, the global influence of the notch was revealed both in the 
deformation of the disk, as indicated by the change in the disk’s maximum radial displacement, and in the 
isolated disk’s center of mass shift. If other factors relating to the single disk system studied here were to 
remain unchanged (e.g., shaft, bearings, etc.), the change in the disk’s center of mass should equivalently 
influence the center of mass shift for the total rotor system. 


Rotordynamic Analysis 

The following rotordynamic analysis was conducted to further understand the dynamic behavior 
of the experimental set-up. In most aviation based turbines, the system is expected to operate 
at some speed above the first critical. The disk spin simulation system being assembled for this continuing 
study consists of a shaft having a diameter of 0.79 in. and a length of 30.75 in., supported by precision 
angular contact ball bearings that are assumed to provide isotropic stiffness. The single disk, as 
seen in figure 1, is mounted at the mid-span of the shaft. The disk’s diametral and polar moments of 
inertia are 63.27 and 126.53 lbm-in 2 , respectively. The weight of the disk is 10.75 lbf. 
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The equations of motion for the rotor system are 


Mq + C(Q.)q + Kq = F(t) 


( 1 ) 


where q denotes the displacement vector; M is the symmetric mass matrix; C(Q) is an asymmetric 
matrix that includes an antisymmetric, speed-dependent gyroscopic matrix; and K is the 
symmetric stiffness matrix. The right hand side of equation (1) represents the force excitation 
vector resulting from the system’s residual imbalance as well as any additional imbalance that 
may result from the disk damage (i.e., notch). 

The system’s natural frequencies are the result of the solution of homogeneous version of 
equation (1), and can be written as 


\ = -w i + Pi 


( 2 ) 


or, in a single degree of freedom form as 


4 



± j<»di 


( 3 ) 


where 0 ) dl = v t is the i th damped natural frequency, and ^ = U ' the corresponding modal damping 

H +v ? 

ratio. The values of co di are directly used for the construction of the whirl speed map (i.e., Campbell 
diagram). The real part of X- t indicates the possible instability of the system. 

Figure 7 is the Campbell diagram for the current experimental system. The intersection 
points of synchronous excitation line IX with the loci of natural frequencies determine the critical speeds. 
It should be noted that the rotor set-up utilized here resembles a Jeffcott rotor [5] with the 
rigid disk located at the center of the flexible shaft, and therefore, the translational and rotational motions 
were statically and dynamically decoupled. Concerning the translational motion, the first two 
natural frequencies of whirl (forward and backward) were independent of the rotor spin speed 
due to the absence of gyroscopic moments. The first synchronous critical speed, i.e., the rotor spin speed 
associated with a resonance condition, was 2610 RPM. During this condition, the maximum shaft 
displacement was at the shaft midpoint as seen in Figure 8. The second forward critical 
speed was at 28153 RPM. This value was beyond the operational speed of the test system which is 
approximately 10000 RPM. 

Figure 9 displays the amplitude and phase of the lx component in the form of Bode plots. 
The first critical at 2610 RPM is indicated by the frequency at which the amplitude peaks and 
a phase inversion occurs. Also, note that for this particular system that the radial displacements 
at the disk location settled to a relatively constant value after the first critical speed was traversed. 
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25 000 






Rotational speed, rpm 


Figure 9. — Bode plot of 1 x component calculated using an assumed 
unbalance of 5 oz. 
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The critical speed map, which plots the critical speeds of the system versus the range of bearing 
stiffness, is presented in figure 10. It is a log -log plot showing the variation of critical speeds for a bearing 
stiffness range of 1x10' to lxlO 7 lb/in. Note that for values above a stiffness of lxlO 7 lb/in., the first 
critical speed did not increase. 

Based on the results of the rotordynamic analysis, the experimental system was correctly assembled 
to operate beyond the first critical speed. Note that the second critical speed was defined by the shaft 
having an “S” shape with zero displacements at the shaft mid-point and the end bearings. Therefore, 
similar to the Jeffcott rotor analysis [5], the system’s synchronous radial displacement at the center of the 
shaft should stay constant beyond the first critical and be a function of the system’s residual imbalance. 
It has been proposed that a cracked disk may have a changing imbalance due to the crack opening as a 
function of rotational speed [2], This is contrary to the assumption of a rigid disk. As a result, the 
amplitude of the synchronous displacement may grow as a function of speed when operating beyond the 
first critical as opposed to being constant. This also will be explored during the experimentation as a 
possible indicator of disk damage. An upcoming analytical study by the authors will address the speed 
dependence of the disk’s eccentricity vector and the influence it may have on the post critical speed 
behavior. 

Based on the results of the rotordynamic analysis, the experimental system was correctly assembled 
to operate beyond the first critical speed. Note that the second critical speed was defined by the shaft 
having an “S” shape with zero displacements at the shaft mid-point and the end bearings. Therefore, 
similar to the Jeffcott rotor analysis [5], the system’s synchronous radial displacement at the center of the 
shaft should stay constant beyond the first critical and be a function of the system’s residual imbalance. 
It has been proposed that a cracked disk may have a changing imbalance due to the crack opening as a 
function of rotational speed [2], This is contrary to the assumption of a rigid disk. As a result, the 
amplitude of the synchronous displacement may grow as a function of speed when operating beyond the 
first critical as opposed to being constant. This also will be explored during the experimentation as a 
possible indicator of disk damage. An upcoming analytical study by the authors will address the speed 
dependence of the disk’s eccentricity vector and the influence it may have on the post critical speed 
behavior. 



Figure 10. — Critical speed map as function of bearing stiffness for 
three lowest critical speeds. 
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Conclusions 


Due to an increased interest concerning the health monitoring of cracked disks in turbine engines, a 
study was required to address the feasibility of utilizing vibration data as a tool for recognizing disk 
damage prior to catastrophic failure. The concept of the approach is based on the fact that the 
development of a disk crack results in a distorted strain field within the component. This may cause a 
measurable deformation in the disk’s geometry and a corresponding change in the system’s center of 
mass. The geometric change and the center of mass shift are indirectly characterized by monitoring the 
amplitude and phase of the first harmonic (i.e., the lx component) of the vibration data. The experimental 
data can be captured by monitoring displacements at either the shaft or blade tips. Because of the limited 
publications concerning the technique and the lack of a basic methodology, the authors undertook the task 
of conducting a modeling campaign as well as preparing the ground work for the experimentation. 

The objectives were met by designing a subscale disk and analyzing the stress states in the 
undamaged state and the notched condition. The disk’s maximum radial displacements and shifts in 
center of mass were defined as functions of notch length and rotational speed. It was seen that the relative 
changes in these parameters were on the order of 0.00014 in. for a 0.962 in. notch in the web area of the 
disk while rotating at 8000 RPM. This value is close to the resolution limits of the displacement sensors 
to be used during the experimentation. The rotordynamic analysis revealed that the experimental system’s 
first critical, assuming a flexible shaft and rigid disk, is 2610 RPM. At this point, it is assumed that the 
experiments will be operated in the range of 0 to 8000 RPM. This allows for plenty of post critical 
operation and at the same time staying within the notched disk’s material limits. 

Upcoming studies from the authors concerning this topic, will include an in depth experimental study 
to verify the above analytical results as well as modifying the Jeffcott rotor analysis to include a speed 
dependent eccentricity vector for the disk. As mentioned above the finite element modeling effort 
indicated relatively small global deformations close to the resolution limits of the existing sensors. The 
subscale experiments, with the inclusion of naturally existing mechanical and electronic noise, will help 
define the true sensitivity of the damage detection methodology. 
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ABSTRACT 

This paper presents the partial results of a dynamic analysis in combination with 
baseline experimental data of a two-spool TF-41, gas turbine aircraft engine. The engine 
was utilized for a government sponsored seeded fault test. By implementing a seeded 
fault in the form of a crack emulating from a notch into the fan stage disk, new sensor 
technologies were assessed in their ability to identify disk cracks in real time. The 
critical speed analysis of the complex engine conducted within this paper was in support 
of the vibration data captured during the full-scale fault tests. Comparisons were 
conducted between the critical speeds predicted by the model and the apparent critical 
speeds obtained from displacement sensors monitoring the fan blade tip clearance in the 
radial direction. Both the analysis and the experimental data indicated multiple critical 
speeds in the domains of operation, although, there were differences in the defined 
critical speeds. It is assumed that discrepancies were introduced by the fact that the 
sensors were monitoring the blades as opposed to the shaft. Therefore, the dynamics of 
the blades may have dominated the actual overall dynamic displacements of the system. 
In addition, disagreements may be due to the modeling assumptions concerning the 
bearing mounts and case rigidity. 

Keywords: Aircraft engine, two- spool rotor, critical speeds, unbalance. 


1. INTRODUCTION 

Recently, the United States set the ambitious goal of reducing the fatal accide nt rate 
for commercial aviation by 80% within ten years [1]. In turn, the National Aeronautics 
and Space Administration (NASA) in collaboration with Federal Aviation 
Administration (FA A) and other federal agencies, universities, as well as airline and 
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aircraft industries, responded by developing the Aviation Safety Program (AvSP). The 
AvSP provides research and technology products needed to help the aerospace industry 
achieve the challenge to improve aviation safety. The Nondestructive Evaluation Group 
at the NASA Glenn Research Center at Lewis Field in Cleveland, Ohio is currently 
addressing the development of propulsion system specific technologies intended to 
detect damage prior to catastrophe under the propulsion health management task. 
Currently, the NDE group is assessing the feasibility of utilizing vibration data for 
detecting cracks in disks. The data is obtained from radial tip clearance measurements 
of blades in turbine engines. The concept of the approach is based on the fact that the 
development of a disk crack results in a distorted strain field within the component. 
This, in turn, causes a minute deformation in the disk’s geometry as well as a possible 
change in the system’s center of mass. 

In order to evaluate emerging realtime sensor technologies for detecting disk cracks 
in turbine engines, the Navy provided an Allison TF-41 turbofan engine. In addition, the 
Navy allowed for the use of a test site at the Naval Air Warfare Center (NAWC) located 
at China Lake, California (see Fig. 1). The TF-41 engine was used by several aircrafts 
during the 1960s and 1970s. It replaced the TF-30 in the A and B model of the Navy A- 
7 aircraft and provided 35% more thrust and a much improved surge margin at high 
power. It is still operational today with the Taiwanese and Greek governments. A cross 
section of the TF-41 engine is shown in Fig. 2. The engine produces about 15,000 
pounds of thrust and consists of two coaxial rotors. For the China Lake tests, a seeded 
fault was implemented in the fan stage disk. The fault was in the form of a notch in the 
web area of the disk. A crack was initiated at the notch tip by running the disk in a spin 
pit facility. Once a crack was initiated and propagated in the spin pit, the disk was 
reassembled in the engine for the full scale tests. 

The objective of the partial rotordynamic analysis conducted in this study was to 
provide a basic understanding of the TF-41’s dynamic behavior. Rotor synchronous 
vibration experienced on this multi- shaft engine results directly from rotor unbalance 
exciting the numerous critical speeds inherent in light-weight, high-speed engine 
structure . Such information was in support of the China Lake experiments that focused 
on using real time vibration data to detect disk cracks. Because of the age of the engine, 
difficulties arose when aiming to acquire original specifications and design data, 
therefore, the need for the modeling presented in this paper. At the time of this writing, 
baseline experimental data was collected at China Lake prior to the implementation of 
the seeded fault. A comparison between the model and the baseline data is provided. 
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Figure 1. TF-41 engine and test stand at China Lake, CA. 



Figure 2. Cross-section of the TF-41 engine. 
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2. DUAL ROTOR -BEARING SYSTEM MODEL 

The requirement of high specific power output for gas turbine aircraft engines has 
resulted in highly flexible rotor designs. These rotors typically operate above several 
critical speeds. The application of rolling element bearings, with low inherent damping, 
makes it difficult to reduce vibrational amplitudes and dynamical loads transmitted to 
the rotor supporting structure. Additionally, the potential dynamic problems are 
amplified by the operation over a wide range of speed and power levels. The typical gas 
turbine engine requires a complex analysis to accurately predict response from rotor 
unbalance throughout the entire engine power range. 

In respect to the TF-41 engine modeled here, be two rotor system consists of 
compressor and turbine stages which are connected by flexible shafting. The inner rotor, 
called the low-pressure rotor (LP rotor), is supported by three rolling-element bearings, 
each of which is mounted to the flexible case structure. The low pressure-intermediate 
pressure compressor (see Fig. 2, [2|) is driven by a two-stage low-pressure turbine. The 
outer rotor, called the high-pressure rotor (HP rotor), is supported principally by three 
rolling -element bearings. There is one inter-shaft differential rolling -element bearing 
connecting the LP and HP rotors. In addition, the LP rotor design does have two spline 
couplings and HP rotor does have one spline joint. 

The schematic of a two-rotor -bearing system is shown in Fig. 3. The two rotors are 
aerodynamically coupled, and have a bilinear low/high speed relationship. 


LP Rotor 



Figure 3. Schematic of computer model of TF41 engine. 


The analysis presented herein includes both flexible rotors, the effect of intershaft 
bearing, and the bearings connecting the rotors to the engine case. The case is 
considered rigid, although casing flexibility may be an important effect in certain 
applications. The analysis includes critical speeds and unbalance response. Since the 
two rotors each have different speeds at a given power setting, a separate critical speed 
analysis is necessary for the frequency range of each rotor. For the given unbalance 
condition, the speed of the exciting (unbalanced) rotor is varied over its entire range. 
The speeds of the excited (balanced) rotor are determined from the provided rotor speed 
relationships [2,3] and determine its gyroscopic effects. Thus, the dynamic analysis 
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considers the spin and whirl speeds of each rotor and defines the LP-excited and HP- 
excited critical speeds. 

The equation of motion for i h rotor can be obtained using Lagrange approach and 
written in the following form: 

M t q t +(C, +cqG )q t + K i q i =F t (t), for i = 1,2 (1) 

where q t is the vector of nodal displacements; M i is the symmetric mass matrix; C, is 
the asymmetric damping matrix including bearings damping; G, is an antisymmetric 
gyroscopic matrix; K i is the asymmetric stiffness matrix including bearings stiffness; 
to, is the rotor speed; and b] is the excitation vector. 

Combining equations (1) for each rotor, the equation of motion for the dual rotor- 
bearing sys tem can be formulated as: 

Mq + ( C + (OjG, + ore, ) q + Kq — F ( t ) (2) 

where M,C and K result from the assembly of matrices M 1 ,M 2 ;C 1 ,C 2 ; and K { ,K 2 , 
respectively. In addition, due to aerodynamic rotor coupling, the speeds are related as 

0)! = aC 0 2 + b (3) 


where a and b are scalars. 


Defining the state vector x T = (xf x 2 ) T = ( q' q 1 ) 7 , the state -space representation of 
homogenous version of Eq. (2) can be written as first order matrix equation 

r o / 

x = Ax , where A — , . , . 

’ -M X K -M~ x ( C + cOjGj + (0 2 G 2 ) 

which can be solved for critical speeds. 


(4) 


3. NUMERICAL SIMULATION RESULTS 

Figure 4 shows the HP rotor exc ited critical speeds (Fig. 4(a)) and the LP rotor 
excited critical speeds (Fig. 4(b)) for the case when all bearings are rigidly mounted to 
the support structure. The summary of critical speed values is presented in Table 1. The 
first two HP-excited critical speeds occur early in the starting transient and the rotors 
will pass rapidly through these speeds conditions. These modes consist primarily of 
pitching and bouncing of the entire engine and pitching of the HP rotor about its spline 
joint. The next two HP-excited critical speeds occur near idle speed and consist of 
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bending of the LP rotor and pitching of the HP turbine about its spline joint. The fifth 
HP-excited critical speed is in the operating speed area and consists of the LP 
compressor pitching about its spline joint and bouncing of two rotors. 



LP Rotor Velocity (RPM) 

(a) 



HP Rotor Velocity (RPM) 

(b) 


Figure 4. Engine critical speeds: (a) HP rotor excited critical speeds, (b) LP rotor excited critical 

speeds. 


The first two LP -excited critical speeds occur near engine idle speed and these 
modes consist of the entire engine pitching and bouncing, and of pitching of the HP 
turbine about its spline joint (in a second mode). The next two LP -excited critical 
speeds occur in the operating range. The mode shapes corresponding to these critical 
speeds consist primarily of bending of the LP rotor and some pitching of the LP and HP 
rotors. 
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HP Rotor Critical Speeds 

LP Rotor Critical Speeds 

(RPM) 

(RPM) 

(431,1500) 

(2100,7319) 

(632,2200) 

(2504, 8211) 

(1920,6680) 

(8255,11776) 

(2246,7818) 

(8710,12060) 

(6377,10524) 



Table 1. Summary of critical speeds arranged in pairs (? v ? 2 ) , where coj is the speed of LP rotor 

and co 2 is the speed of HP rotor. 


4. PRELIMINARY EXPERIMENTAL RESULTS 


Figure 5 illustrates the magnitude and phase of vibration responses from three 
separate sensors; one eddy current sensor and two capacitive sensors. Each sensor is 
aligned to monitor the radial blade tip clearance of the TF-41 fan stage. The data was 
captured during full scale engine tests prior to the implementation of the flaw. To 
further accentuate the critical speed response of the engine, a 1.06 ounce unbalance at a 
radius of 5.6 inches was added to the 32 inch fan disk [5]. The procedure for calculating 
unbalance amplitude and phase involved a least mean squares regression of a 
1 /revolution sinusoid to the interblade spacing. The following relationship [6] was 
developed from a simple kinematic analysis of a rotor with evenly- spaced blades 
spinning about a center which is not its geometric center: 


where: A R 
R 


A0 

0 

avg 


AR = 


A0 R 



avg 


is the unbalance amplitude 
is the tip radius 

is the complex amplitude of the 1/rev sinusoidal fit to interblade spacing 
is the average interblade spacing 


The experimental data of Hg. 5 suggests multiple critical speeds with the most 
obvious ones at approximately 6,300 and 8,000 rpm. The values for the experimentally 
obtained critical speeds do not appear to coincide with all the numerically predicted 
values as seen in Table 1. 


The discrepancies between the numerical predictions and the experiments may be 
due to multiple sources. First, focusing on the experiments, it is assumed that a portion 
of the discrepancies can be attributed to the fact that the sensors were monitoring the 
blades as opposed to the shaft. As a result, the motions of the blades may have 
overridden the actual dynamic displacements of the system. Another source of the 
disagreement may be due to assumptions made when creating the numerical model. The 
bearing mounts were assumed to be rigid. The experimental data may suggest dampers 
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at some of the bearing mounts. Lastly, due to a lack of information, the engine case was 
modeled as a rigid structure. 



RPM 

Figure 5. Unbalance magnitude and phase calculated from interblade spacing. (30gms of 
unbalance mass added at a 5.6” radius for this run) [4]. 


5. CONCLUSIONS 

In summary, this study provided a numerical dynamic analysis and compared the 
results with the available baseline experimental data concerning a two-spool TF-41, gas 
turbine aircraft engine. The engine was utilized for a government sponsored seecbd fault 
test. The overall goal of the full scale testing was to assess new sensor technologies for 
detecting cracks in an actual fan disk during engine operation. The critical speed 
analysis of the complex engine conducted within this paper was in support of the 
vibration data captured during the full-scale fault tests. 

Both the numerical model and the experimental data indicated multiple critical 
speeds in the region of operation. When comparing the modeling results with the data, it 
was apparent that there were differences in the defined critical speeds. It is assumed that 
these discrepancies were introduced by the fact that the sensors were monitoring the 
blades as opposed to the shaft. The motions of the blades may have overridden the 
actual dynamic dsplacements of the system. In addition, discrepancies may be due to 
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the modeling assumptions concerning the rigidity of the bearing mounts as well as the 
engine case. 
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ABSTRACT 

The transient vibration response of a cracked flexible rotor passing through its critical speed is analyzed for crack 
detection and monitoring. The effects of different factors such as various crack depths, acceleration, damping, torque, 
unbalance eccentricity, and rotor weight on the rotor vibrational response are studied. The breathing types of cracks are 
analyzed using simple hinge model in a case of shallow cracks, and the cosine function is employed in the case of deep 
cracks. The developed strategy enables the analysis of cracked rotor vibrational response with and without weight 
dominance, taking into account also the nonsynchronous rotor whirl. In addition, the local cross-flexibility for deep 
cracks is taken into account. Lastly, the effect of the crack depth on “stalling” of the rotor due to the limited driving 
torque is investigated. 

Keywords: rotor, crack, acceleration, vibrational response, detection, monitoring. 

1. INTRODUCTION 

In many cases the rotors of modern machines are rapidly accelerated from rest to operating speed to reduce the excessive 
vibrations at the critical speeds. Also, the vibration monitoring during startup or shutdown has been receiving growing 
attention, especially for machines such as aircraft engines, which are subjected to frequent starts and stops as well as 
high speeds and acceleration rates [1]. It has been recognized that the presence of angular acceleration strongly affects 
the rotor’s maximum response to unbalance and the speed at which it occurs. 

The topic of transient cracked rotor response has been treated by number of published works [2-4]. They have been 
focused on studying the dynamic behavior of a rotor with the so-called breathing crack while passing through a critical 
speed at constant angular acceleration or deceleration. Sekhar [2-3] investigated the transient vibration response of a 
cracked rotor passing through the critical speed utilizing a simple hinge model for small cracks. He followed an 
assumption that the vibrations remain small in comparison to the rotor’s sag. If a cracked shaft rotates under the load of 
its weight, then the crack will open and close once per revolution, in a case of synchronous whirl. Sawicki et al. [4] 
studied accelerating cracked rotor response using the angle between the crack centerline and the rotor whirl vector to 
determine the closing and opening of the crack. This allowed for studying the rotor dynamic response with and without 
the rotor weight dominance by taking into account nonsynchronous whirl. 

A few authors [5-7] have studied the dynamic response of uncracked rotors subjected to a constant driving torque. 
Markert et al. [5] determined the maximum deflection of a Jeffcot rotor during the acceleration through the critical speed 
or during stalling in resonance as function of the three nondimensional system parameters, i.e., damping ratio, 
eccentricity ratio and driving torque. Later, Gasch et al. [6] applied these strategies to a practical flexible rotor with a 
continuous mass distribution passing through its critical speeds while subjected to a constant driving torque. Genta and 
Delprete [7] developed a mathematical finite element-based model to study the nonlinear behavior of complex 
anisotropic rotors with non-constant angular speed. 
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There is a wealth of published research results on cracked rotors. Bendy [8] presented a rotor crack detection procedure 
based on the start-up and shutdown vibration monitoring. Wauer [9] provided a comprehensive literature review 
concerning the state-of-the-art vibration methodologies related to cracked structures. Gasch [10] provided a 
comprehensive investigation of the stability behavior of a cracked Jeffcott rotor and the forced vibration due to 
unbalance and crack at constant rotational speed. Penny and Friswell [11] demonstrated the influence of the crack model 
on the response of the Jeffcott rotor. 

In this paper the effect of a breathing crack on the unbalance response of a accelerating rotor subjected to a constant 
driving torque is studied. The equations for the rotor’s lateral motion are supplemented by an additional equation relating 
the angular displacement and the driving torque. Torsional elasticity of the shaft is neglected due to the assumption that 
the torque is applied to a torsionally stiff rotor. In addition, the behaviors of a stalled rotor’s angular velocity as well as 
the rotor’s maximum amplitude of unbalance are presented as functions of the crack characteristics. The consequences of 
torsionally induced slant cracks are not considered in this study. 


2. CRACK MODELS AND STIFFNESS MATRIX 

The theoretical model, called the Jeffcott rotor, employs a flexible rotor composed of a centrally located, 
unbalanced disk attached to a massless elastic shaft. In turn, the shaft is mounted symmetrically on rigid bearings (see 
Fig. 1(a)). For this study, the shaft has a transverse crack running across its cross-section located near the disk. The 
stiffness of the uncracked rotor system is symmetric (isotropic) and the damping due to the air resistance effect is 
assumed to be viscous. The rotor is driven by a constant external drive torque. The angle between the crack centerline 
and the line connecting the bearings and shaft center (Fig. 1(b)), y/ = <J> - arctan (j/z) , is used to determine the closing 

and opening of the crack. At any instant of time, the c-axis remains perpendicular to the face of the crack. This allows 
the body-fixed rotating coordinate frame ( C, rj, c) to spin with the same velocity as the rotor. The weight dominance, 
assumed in almost all previous analyses of horizontal cracked rotors, is not required. Furthermore, the influence of the 
whirl speed on the closing and opening of the crack is included. This nonlinear system with time-varying stiffness 
coefficients is studied numerically, with particular focus on the effect of different crack depths on the rotor stalling. 




Figure 1 . Model of Jeffcott rotor with crack and constant driving torque (a), crack section in inertial and rotating coordinates (b). 


The stiffness matrix for a Jeffcott rotor with a cracked shaft in rotating coordinates can be written as: 
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where the first matrix refers to the stiffness of the uncracked shaft, and the second defines the changes in stiffness Ak 
and A k t] in c and // directions, respectively. The function /((//) is a crack steering function which depends on the 
angular position of the crack y/ and the selected crack model. 

The simplest representation of a crack is defined as the hinge model. In this representation, the crack is assumed to 
change from its closed to open state suddenly as the shaft rotates. The steering function for this model is defined as 


or, using Fourier expansion is 
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While the hinge model might be proper representation for very small cracks, Mayes and Davies [12-13] proposed a 
model with a smooth transition between the opening and closing of the crack, which is more appropriate for larger 
cracks. In this case the crack steering function, or the Mayes modified function, takes the following form: 

l + cos(^) 
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(3) 


Next, the stiffness matrix for a Jeffcott rotor with a cracked shaft in inertial coordinates, K y is given by 

K 7 = TK^T -1 

where the transformation matrix T is 


T = 


cosd> -sin® 
sin® cos® 


Thus 


K, = TK fi T 1 = 



' K O'' 

I f(¥)K 

f Ak\ + Ak 2 cos 2® A k 1 sin 2® ' 


,0 

1 2 

v A k^ sin 2® A/q - Ak 1 cos 2® y 


(4) 


(5) 


(6) 


A k, + A k„ Ak, - A k„ 

where A h = — 2 - and A k 2 = — - . 
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3. EQUATIONS OF MOTION OF CRACKED ROTOR UNDER CONSTANT DRIVING TORQUE 


The nonlinear coupled equations of motion for the accelerating Jeffcott rotor with a cracked shaft, subjected to constant 
driving torque, T a , unbalance force, and gravitational force due to its weight, can be written in the inertial coordinate 
frame as: 



(7.1) 
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, are the rotor mass and damping matrices, respectively; q = (z yf 


The matrices, M = 


the vector of the disk’s displacements; and d? is the rotor spin angle (see Fig. 1(b)). 


is 


Using Eq. (7.1) one can obtain the following equation relating torque and the angular displacement: 

J p <b = T a + sMg sin 6 + sK ( y cos 6 - z sin ) + sC ( y cos 9 - z. sin 0) + x 

2 (7.2) 
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where # = d?-/? , y/ = Q>- arctan^— j, # = d>, # = <!>, and J p =Mr 2 . 

Normalizing displacements with respect to unbalance eccentricity, and introducing nondimensional time, damping, 
and torque, the equations (7.1) and (7.2) take the following nondimensional form: 
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The following definitions for nondimensional variables were employed in Eqs. (8.1) and (8.2): 
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4. RESULTS AND DISCUSSION 


Five nondimensional parameters were used in the numerical study of the accelerated cracked rotor, i.e., T, (ff/r) , A K , 

g 

C , and — — . The following results were generated for the case where unbalance is in the direction of the crack, i.e., 


SCO,, 


Ak r 


P = 0 . Also, the cross-stiffness for deep cracks was represented as Ak = — — [13]. 
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Figure 2. Normalized vibration amplitude (a) and normalized spin speed (b) as a function of time; 
7 = 0.01 ,g = 0.02, (s/rf = 0.85 x 10~ 3 , andgy^®,?) = 5 . 


Figure 2 illustrates the behavior of the rotor passing through the critical speed for the uncracked and cracked rotor 
under the constant driving torque, T = 0.01 . The normalized rotor vibration amplitude and normalized spin speed are 
shown as a function of nondimensional time. Note that the appearance of a crack made the zone of critical speed wider. 
Next, the slight increase of the depth of the crack, here by 0.001, caused the rotor to achieve a stalled condition, i.e., it 
failed to accelerate beyond the critical speed. In other words, all the power delivered by the given torque was dissipated 
by nonrotating damping [7] and the rotor acceleration was no longer possible. 
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Figure 3. Normalized shaft deflection and spin speed as a function of time for g = 0.02 (a), g = 0.04 (b); 
AK=0.1, r = 0.01, (f/r) 2 =0.85xl0~ 3 , g/(«u 2 ) = 5. 


The effect of damping is illustrated in Fig. 3. The increased damping suppressed the stall mechanism and allowed 
the rotor to traverse through the critical speed. 

Figure 4 illustrates that the appearance of crack caused 1/3 and 'A subharmonic peaks to show up and also 
significantly increased the amplitude of vibration of the “stalled” rotor. In this case the large vibration amplitudes of the 
stalled cracked rotor might well exceed the rotor static deflection and therefore violate the common weight dominance 
assumption made in the past studies concerning breathing cracks. 
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Figure 4. Time history of normalized vibration amplitude for uncracked rotor (a) and with crack A K = 0.25 (b) at the stalled 

condition; g = 0.03, T = 0.001, (£•/ r) 2 =0.001, gy/|£<w 2 j = 5 . 


For the rotor in stall and under the constant driving torque, the normalized spin speed and vibration amplitude as a 
function of crack depth is presented in Fig. 5. It can be seen that while the rotor vibration amplitude increased 
parabolically (Fig. 5(a)), the stall speed showed a near linear dependence on the depth of the crack (Fig. 5(b)). 
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Figure 5. Normalized spin speed (a) and maximum vibration amplitude (b) as a function of crack depth at the rotor stalling; 

T = 0.01, (f/r) 2 =0.85xlCT 3 , g/(^j = 5, £ = 0.02. 


The effect of the applied constant torque and various crack depths on the maximum rotor resonant vibration 
amplitude and the spin speed for various crack depths is shown in Fig. 6. The peak of each curve in Fig. 6(a) denotes the 
threshold value of the constant driving torque below which the rotor is locked in the resonance and above which the 
rotor is capable of surpassing the resonance. 
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Figure 6. Maximum normalized vibration amplitude (a) and the corresponding maximum normalized spin speed (b) of the rotor as a 


function of constant torque for various crack depths; £ = 0.03, (e/r) - = 0.001 , gj(s(0 2 j 


= 5 . 


For the range of torques below the threshold value (i.e., in the stall zone), the rotor’s maximum vibration amplitude 
increased with increasing torque (Fig. 6(a)). Interestingly, the stall speed remains almost constant (Fig. 6(b)), indicating 
a lack of rotor acceleration. Also, while in the stall zone, the maximum vibration amplitude grew with increased depth of 
the crack, and the corresponding rotor stall speed decreased. By increasing the drive torque, the cracked rotor is allowed 


to accelerate through resonance, reduce its maximum resonance vibration amplitude, and shift the resonance peak to 
higher spin speeds. 


5. CONCLUSIONS 

The following conclusions can be drawn based on the results presented in this paper: 

1 . The developed model enables the study of the cracked rotor dynamic response with and without weight dominance 
in the presence of nonsynchronous whirl. The model includes small and deep cracks, as well as cross stiffness effect 
(for deep cracks). 

2. The stalling effect of the rotor is crack sensitive, i.e., even a minute increase in crack depth can cause the rotor to 
stall. This is accompanied by the presence of subcritical (1/3, V-i) response peaks and a large increase in the 
fundamental vibration response. 

3. For a rotor in stall, the maximum vibration amplitude increases for increasing value of torque while the 
corresponding stall spin speed remains almost constant. 

4. In addition, the maximum vibration amplitude during stall grows parabollically with increased crack depth. The 
corresponding rotor stall speed decreases almost linearly with crack depth. 
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external damping coefficient 
crack steering function 
moment of inertia of the disk 

stiffness of uncracked shaft 
mass of the disk 
radius of gyration 
time; r = co n t 

driving torque; T =T a /j p col 
inertial coordinates; Z = z/s, Y = y/e 
acceleration ratio 

angle between crack and unbalance eccentricity 

the largest stiffness change ratio in q -direction caused by crack 

the largest stiffness change ratio in rj -direction caused by crack 

stiffness change ratio (= A kg IK j 

eccentricity of the disk 
external damping ratio 

angular position of the eccentricity s in the inertial system 
body fixed rotating coordinates, £ is in the crack direction 
undamped critical speed of uncracked rotor 
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ABSTRACT 

Rotor crack problems present a significant safety and loss 
hazard in nearly every application of modern turbomachinery, 
particularly in the power generation industry. However, early 
crack detection is not easily achieved during the operation of 
machinery. The difficulty is based on the fact that a crack 
produces an undetectable change in the overall structural 
response. This paper analyzes the coupling of torsional and 
lateral vibrations for an unbalanced cracked rotor. The rotor 
equations of motion for a system with cracked shaft, obtained 
using Lagrangian dynamics, show coupling and nonlinear 
interaction between the torsional and lateral vibrations. To 
investigate the effect of a transverse surface crack on the 
dynamic rotor response the breathing crack model was 
employed. By applying an external torsional excitation together 
with the excitation due to unbalance, signature responses were 
observed in the rotor vibration spectrum at sum and difference 
frequencies. These signature responses were due to the 
nonlinear effect of the crack. The observed phenomena, 
analytically defined here, offers a new methodology concerning 
crack detection and prognosis in rotors. 

INTRODUCTION 

In many cases of rotating machinery the rotors have the 
lowest torsional natural frequency higher than the second or 
third bending natural frequency. The lateral rotor vibrations are 
always an important consideration for virtually all types of 
rotating machinery. Conversely, the torsional rotor vibration is 
often not an important consideration for many machinery types, 
especially machines with single uncoupled rotors. In addition, 
the measurement of torsional vibrations, even those having 
amplitudes close to the destructive level, is rather difficult task, 
since complex and specialized transducers need to be 
employed. Consequently, the attention paid to the monitoring 


of torsional vibration, as well as to the dynamic coupling 
between the torsional and lateral rotor vibration characteristics 
is typically discounted. 

A study by Tondl [1] was one of the first to investigate the 
effect of the coupled lateral and torsional vibrations on 
turbogenerator rotor stability. He concluded that the combined 
effect of torsional stiffness and bending stiffness results in 
speed intervals where the rotor vibrations due to residual 
unbalance become unstable. Loewy and Piamlli [2] also 
analyzed the problem of coupled vibrations, but due to 
incorrect derivation of the rotor equations of motion they 
arrived at the wrong conclusions. 

The online monitoring of machine dynamic behavior 
cannot be overestimated. Any new advancement in detection 
and diagnosis of critical malfunctions in rotating machinery can 
be extremely beneficial to industry [3]. This especially applies 
to fatigue cracks within the shaft, which present a potential 
source for catastrophic failures of rotating machinery. Most 
reported studies are focused on two crack signatures, i.e., twice 
the running frequency component (2X) and the subharmonic 
component at approximately half of the shaft critical speed. 

Papadopoulos and Dimarogonas [4] used a non-rotating 
cracked Timoshenko shaft to demonstrate the existence of an 
apparent coupling of torsional and bending vibration. They 
modeled crack using the local flexibility matrix, and then 
proceeded to study the vibration spectra in the presence of 
harmonic excitations. Other researchers [5-9] also dealt with 
the problem of coupled vibration. Recently, Muszynska et al. 
[5] analytically and experimentally analyzed lateral/torsional 
coupling mechanisms resulting from combinations of the 
unbalance, shaft stiffness asymmetry, and radial sideload. 
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Muszynska experimentally observed torsional resonance at 
speeds equal to 1/8, 1/6, 1/4, and 1/2 of the lowest torsional 
natural frequency. Bently et al. [6] continued this study with 
special attention paid to the analysis of a “snapping” action that 
occurs when during rotation the rotor experiences a peak of 
torsional acceleration. 

The topic of cracked rotor vibrations has been analyzed in 
a number of published works [9-13]. They have been focused 
on the study of dynamic behavior of rotors with the so-called 
breathing type of crack during the passage through a critical 
speed at constant angular acceleration or deceleration. For 
example, Sawicki et al. [10] studied the accelerating cracked 
rotor response using the angle between the crack centerline and 
the rotor whirl vector, to determine the closing and opening of 
the crack. This allows one to study the rotor dynamic response 
with or without the rotor weight dominance by taking into 
account nonsynchronous whirl. Gasch [12] provided a 
comprehensive investigation of the stability behavior of a 
cracked Jeffcott rotor. The cracked shaft analysis of Gasch 
assumed a constant rotational speed and a forced vibration 
resulting from the system’s residual unbalance. 

In this paper, an analytical solution of a Jeffcott rotor with 
a breathing crack is developed in order to show coupling 
between the torsional and lateral vibrations. The equations for 
the rotor’s lateral motion are supplemented by the third 
equation for torsional vibration. Though the assumed model for 
the cracked rotor is based on an idealized Jeffcott rotor, the 
generated results provide a new approach for understanding the 
complex vibration response of cracked rotors. 

NOMENCLATURE 

C L Rotor lateral damping 

C T Rotor torsional damping 

F , , F External forces 

f(ys) Crack steering function 

J Moment of inertia of the disk 

K y Shaft stiffness matrix in inertial coordinates 

Ky; Shaft stiffness matrix in rotating coordinates 

K Lateral stiffness of uncracked shaft 
K t Shaft torsional stiffness 

M Mass of the disk 

z, y Inertial coordinates 

A kc Stiffness change ratio in q -direction due to crack 

Ak t) Stiffness change ratio in rj -direction due to crack 

A K Stiffness change ratio A kg IK j 

T Transformation matrix 

T e External torque excitation 

t Instantaneous time 

1 3 Angle between crack and unbalance vector 

e Eccentricity of the disk 

Rotor-fixed rotating coordinates 
, C T Lateral and torsional damping ratio 
0 Angular position of the eccentricity vector 


<t> Angle between If and z axes 

y/(t) Torsional angle 

co Rotor angular speed 

w e Torsional excitation frequency 

CRACKED ROTOR MATHEMATICAL MODEL 

The theoretical model, called the Jeffcott rotor, employs a 
flexible rotor composed of a centrally located unbalanced disk 
attached to a massless elastic shaft which, in turn, is mounted 
symmetrically on rigid bearings. Here, the shaft has a 
transverse crack running across its section and located adjacent 
to the disk. The sectioned view of the cracked rotor, in both 
inertial and rotating coordinates, is shown in Fig. 1. The 
stiffness of the uncracked rotor system is symmetric (isotropic), 
and the damping due to the air resistance effect is assumed to 
be viscous. The angle between the crack centerline and the line 
connecting the bearings and shaft center, y/ = <J> - arctan (j/z) , 
(see Fig. 1) is used to determine the closing and opening of the 
crack. At any instant of time, the c-axis remains perpendicular 
to the face of the crack, causing the body-fixed rotating 
coordinate frame ( C >]■ C) to rotate with the same velocity as the 
rotor. 


in 



Figure 1. Rotor crack section in inertial and rotating 
coordinates. 


The angular position of the unbalance vector can be 
expressed as 0(t) = cot + y/{t) + ft , where co is a constant spin 
speed of the shaft, y/(t) is a torsional angle, and ft is the fixed 
angle between the unbalance vector and the centerline of the 
transverse shaft surface crack. The kinetic and potential energy 
for the rotor system subjected to lateral and torsional vibrations 
can be expressed as follows: 

T = -^-M jz 2 + y 2 + e 2 0 2 +2s6\-zsmd + ycos6^+^J p 6 2 (1) 
U = \{z j}Kyj^J + ii^ 2 (2) 

where 
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is the shaft stiffness matrix, with the off-diagonal elements 
corresponding to asymmetry induced either by a crack or the 
geometry of the system. Utilizing the Lagrange approach, the 
nonlinear, coupled equations for transverse and torsional 
motion for the rotor system take the following form: 

Mz + C L z + K zz z + K. v y = F z + M£- ( # 2 cos <9 + <9 sin 

V (4) 

My + C L y + K yz z + K yy y = F y + Mtt ( <9 2 sin d-d cos 6 »j 

J p y/ + C T y/ + sC L (isin#- ycosd)- s(F z sin 6-F y cos#) + 
s \[ K zz z + K zyy ) sin 6 - ( K y: z + K^y} cos <?] + = T e (t) 


K ; = TK r T = 


X V 

K „J 


'K 0 ) 

f(¥)K I 

'A k[ + Ak, cos 2® 

Ak, sin 2® ) 

,0 K) 

2 1 

Ak, sin 2® 

Ak, - Ak, cos 2 ® J 


( 10 ) 


where the transformation matrix T is 
^cos® -sin® 


T = 


sin® cos® 


and 


Ak, + Ak Ak, - Ak 

At = n - and Ak, = ■ 1 " 

1 T7 - L 


K 


K 


( 11 ) 


( 12 ) 


Uncracked (symmetric) shaft 

For the uncracked rotor the following conditions related to 
potential energy and stiffness components are satisfied: 

J^ = K tV, K zz = K yy ^ K and K zy = K yz = 0 ( 6 ) 

Cracked shaft 

The stiffness matrix for a Jeffcott rotor with a cracked shaft 
in rotating coordinates can be written as: 



(K 0 ^ 

-f(v) 

( Ak, 

0 2 

K*=| 

v 0 K } 

0 

V 



( 7 ) 


where the first matrix refers to the stiffness of the uncracked 
shaft, and the second defines the changes in stiffness Ak, and 

A k in £and 77 directions, respectively. The function /((//) is 

a crack steering function which depends on the angular position 
of the crack, y/ , and the selected crack model. 


The simplest crack model is the hinge model, where the 
crack is assumed to change from its closed to open state 
suddenly as the shaft rotates. The steering function for this 
model is defined as: 



0 for £ < 0 

1 for £ > 0 


( 8 ) 


The cross-stiffness for deep cracks has been accounted for 

Akc 7 Akc 5 Akc 

as A k„ = — — , which yields At = — and Ak, = — . 

n 6 1 6 K “ 6 K 


The derivative of the potential energy term in Eq. (4) 

dy/ 

takes the following form: 


du _ K 77 (®) 

dy/ 4 dy/ 

|j Akj + Ak 2 cos 2® ) z 1 + 2zyAk 2 sin 2® + ( Ak, - A k 2 cos 2® ) y 2 J 
+ ,/(®)AiAk 2 |~^ 2 s j n 20 - 2zv cos 2® - y 2 sin 2® j + K T y/ 

( 13 ) 


NUMERICAL SIMULATIONS 

To investigate coupling of the rotor bending and torsional 
vibrations induced by lateral excitation (unbalance) and/or a 
rotor’s crack, a number of numerical simulations were 
conducted. The realistic physical parameters, listed in Table 1, 
were selected for the model, in order to highlight the relevance 
of the study to real life applications. Runge-Kutta integration 
was applied to obtain numerical solutions of the nonlinear 
coupled equations of motion (4)-(5). The integration time 
sampling frequency was selected as 200kHz. 

Table 1: Numerical Model Physical Parameters. 


While the hinge model might be an appropriate 
representation for very small cracks, Mayes and Davies [12-13] 
proposed a model with a smooth transition between the opening 
and closing of the crack that is more adequate for larger cracks. 
In this case the crack steering function, or the Mayes modified 
function, takes the following form: 

l + COS((77) 


f(w ) = - 


( 9 ) 


The stiffness matrix for a Jeffcott rotor with a cracked shaft 
in inertial coordinates, K ; is derived as 



Physical parameter 

Value 

Units 

M 

Disk mass 

20 

kg 

J 

Disk polar moment of inertia 

0.3 

kg m 2 

K 

Shaft stiffness 

8.82xl0 5 

N/m 

s 

Unbalance eccentricity 

1.5xl0" 3 

m 

P 

Unbalance phase angle 

0.0 

degrees 

4l 

Lateral damping ratio 

0.01 


4t 

Torsional damping ratio 

0.001 


co L 

Lateral natural frequency 

210 

rad/sec 

w T 

Torsional natural frequency 

630 

rad/sec 


3 


Copyright © 2004 by ASME 



Figure 2 illustrates the vibrational spectra for the 
uncracked rotor running at the constant speed <y=100 rad/sec, 
i.e., below the bending critical speed, and subjected only to 
unbalance and gravitational forces. The nonlinear dynamic 
coupling between the lateral and torsional modes can be seen 
clearly in the torsional response spectrum (Fig. 2(a)), where the 
sum and difference of lateral natural frequency ( w L ) and the 
running speed (co), co L ±co, are centered with respect to lateral 
natural frequency. In addition, the torsional natural frequency 
component is present. It can be shown that when the running 
speed is higher than the bending critical speed, the two 
frequency components co±co L show up and are located 
symmetrically with respect to running speed. 



(a) Torsional vibration spectrum. 

xIO ' 3 



(b) Lateral vibration spectrum. 

Figure 2. Vibrational spectra due to unbalance for uncracked 
rotor at speed ry=100 rad/sec. 


The corresponding lateral vibration spectrum (Fig. 2(b)) 
shows primarily a strong synchronous frequency component 
and small amplitude peak corresponding to the bending critical 
speed. 


Next, a harmonic torsional excitation, T e {t) = 2sin(210t) , 
is applied to the rotor, in order to investigate its effect on lateral 
vibration in a context of coupling between torsional and lateral 
vibrations. This torsional excitation is in addition to the 
unbalance and gravitational forces. From Fig. 3(a) it can be 
seen that the frequency content and the amplitude of torsional 
vibrations show negligible change (as compared with Fig. 2(a)), 
except the appearance of the torsional excitation frequency of 
210 rad/sec. The torsional excitation has no affect concerning 
the lateral vibration spectrum, as seen in Fig. 3(b). 

It can be concluded that for the case of an uncracked 
Jeffcott rotor (supplemented with the third equation of torsional 
vibration), the existing unbalance excitation in the lateral 
direction generates a torsional response, while the effect of 



(a) Torsional vibration spectrum. 



(b) Lateral vibration spectrum. 

Figure 3. Vibrational spectra due to unbalance and torsional 
excitation T e (t) = 2sin(210t) for uncracked rotor at 
speed (o= 100 rad/sec. 

applied torsional excitation on the lateral vibrations is 
negligible. 
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The next step involves the case of a cracked rotor subjected 
to forces related to unbalance and gravity (i.e., no torsional 
excitation). The assumed cracked depth is A K = 0.4 and the 
rotor speed is <y= 21 rad/sec. The torsional and lateral frequency 
vibration spectra are shown in Fig. 4. It can be noted that both 
the torsional and lateral spectrum indicate IX, 2X, and 3X 
frequency components. This is due to the nonlinear mechanism 
of the crack. The coupling between the torsional and lateral 
vibrations is manifested by the presence of synchronous 
frequency and its higher harmonics in the torsional spectrum 
(see Fig. 4(a)) in the absence of torsional excitation. Note that 
for the uncracked rotor, there were no synchronous components 
in the torsional vibration spectrum (see Fig. 2a). 


frequency, are generated by the coupling between the torsional 
and lateral vibrations, induced by the presence of the crack in 
the rotor. 

Finally, a torsional excitation of T e (t) = 200sin(210t) was 
applied to the case of the cracked rotor, again, in addition to the 
unbalance and gravity forces. These results are shown in Fig. 5. 
Figure 5(a) indicates that the rotational frequency component 
and its higher harmonics are suppressed by the large amplitude 
of the harmonic torsional excitation. The only visible responses 
are at the torsional excitation frequency and the torsional 
natural frequency. A different scenario is presented by the 
rotor’s lateral vibration spectrum, as shown in Fig. 5(b). In 
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(a) Torsional vibration spectrum. 



(b) Lateral vibration spectrum. 

Figure 4. Vibrational spectra due to unbalance for cracked 
rotor at speed aj= 21 rad/sec. 


The detailed analysis of the torsional vibration spectrum in 
Fig. 4(a) reveals, in addition to torsional resonance frequency 
co T = 630rad/sec., the presence of the sum and difference 
frequencies {a> L ±a>/ 2), ( a> L ±3a>/2 ), ( a> L ±5a>/2 ), and 

{a> L ±la>/2), which are located around the bending natural 
frequency ( co L ). These frequencies, being the result of 
interaction between the lateral natural frequency and rotational 


(a) Torsional vibration spectrum. 


xIO ' 4 



(b) Lateral vibration spectrum. 

Figure 5. Vibrational spectra due to unbalance and torsional 
excitation T e (t) = 200sin(210t) for cracked rotor at 
speed co= 21 rad/sec. 

addition to such crack “indicators” as the 2X and 3X 
components, the torsional excitation frequency ® c =210 
rad/sec. is seen along with the side frequencies w e ± w , 
co e ±2co, and co e ± 3® , which are located around w e . These 
frequencies are generated as a result of the non-linear 
mechanism of the breathing crack under the influence of 
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imbalance, gravity, and harmonic torsional excitation. Again, 
the differences concerning the torsional and lateral vibration 
spectrums between the uncracked and cracked cases are quite 
obvious as discussed above. 

CONCLUSIONS 

An analytical solution was developed for the Jeffcott rotor 
taking into account coupling between the torsional and bending 
vibrations. This was achieved by the inclusion of the third 
equation for torsional vibration. The modeling was further 
enhanced by considering the nonlinearities associated with a 
breathing crack located on the shaft adjacent to the disk. 
Numerical simulations of the uncracked and cracked cases were 
carried out by applying an unbalance excitation, gravity forces, 
and a torsional excitation. 

When comparing the torsional and lateral frequency 
spectrums for the uncracked and cracked cases, very clear 
distinctions were noticed. In each case (uncracked and 
cracked), the system was excited by the residual unbalance and 
gravitational forces, and then by the addition of a torsional 
excitation (i.e., 4 cases total). The sum and difference 
frequencies have been observed in a response of cracked rotor 
due to unbalance alone and due to unbalance and external 
torsional excitation. In the case of unbalance excitation only, 
the torsional spectrum exhibits the sum and difference 
frequencies located symmetrically around the bending critical 
speed. In a case of applied torsional excitation, in addition to 
the present unbalance, the lateral vibration spectrum shows sum 
and difference frequencies located symmetrically around the 
bending critical speed. 

The uniqueness of the observed crack signatures, captured 
by the analytical approach, demonstrates an improvement in 
early diagnosis of rotor cracks. 
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ABSTRACT 

The nonlinear model of the cracked Jeffcott rotor is investigated, with the particular focus on study of rotor’s vibrational 
response using tools of nonlinear dynamics. The considered model accounts for nonlinear behavior of the crack and 
coupling between lateral and torsional modes of vibrations. Load torque is applied to the rotor which is laterally loaded 
with a constant radial force (gravity force) and unbalance excitation. The co-existence of frequencies of lateral modes in 
the frequency spectra of torsional mode are characteristics of the coupling response of lateral and torsional vibrations. 
When only the lateral excitations are applied, vibration amplitude bifurcation plot with the shaft speed as a control 
parameter, demonstrates some speed ranges for which vibrations of the rotor dramatically increase. Furthermore, the 
torsional response amplitude at the same speed ranges also increases and chaotic behavior can be observed due to the 
lateral excitations. These phenomena cannot be observed for pure lateral vibration response with the torsionally rigid 
rotor assumption. 

Keywords: Jeffcott rotor, crack, nonlinear dynamics, vibration signatures, bifurcation diagram, detection. 

INTRODUCTION 

The tendency to higher speeds in turbomachinery results in design of more flexible shafts, which spin at speeds above 
several of their natural frequencies. Due to stress concentration and high spin speed, rotor dynamic system is more 
prone to cracks which propagate faster due to low-cycle fatigue loading. In many rotating machinery, the lateral natural 
frequencies are far much lower than the lowest torsional natural frequency. That’s why most of the condition 
monitoring methods is based on lateral vibrations. Without proper transducer, it’s very difficult to detect torsional 
vibrations even when their amplitudes are close to dangerous level. However, in order to capture the whole picture of 
the cracked dynamic system, and develop a high efficiency on on-line condition monitoring method, the lateral and 
torsional coupling must be taken into consideration, because torsional excitation generate not only torsional responses, 
but also cause lateral vibrations due to the mechanism of torsional and lateral coupling. 

Finite Element Method is basically a linear approach, which ignores the lateral and torsional coupling mechanism of the 
uncracked rotor. It is only due to the crack, which introduces nonlinearity and produces coupling terms of cracked rotor 
in the FEM-based model. But, the coupling does exist in an uncracked rotor dynamic system. Tondl [1], Cohen and 
Porat [2] and Bernasconi [3] concluded that typical lateral excitations, such as unbalance, may result in both lateral and 
torsional responses of uncracked rotor. Furthermore, Muszynska et al. [4] and Bendy et al. [5] discuss rotor coupled 
lateral and torsional vibrations due to unbalance, as well as due to shaft anisotropy and radial constant preload force. 
Their experimental results exhibited the existence of significant torsional vibrations, due to coupling with the lateral 
modes. Sawicki et al. [6,7] studied transient response of the cracked rotor, including the stalling effect under the 
constant driving torque. 

EQUATION OF MOTION OF A CRACKED JEFFCOT ROTOR 

A cracked Jeffcott rotor with lateral and angular (torsional) degrees of freedom is considered. The rotor is subjected to a 
constant radial load (gravity force), rotating unbalance force, and the externally applied torque. The coupled torsional- 
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flexural vibrations are studied for a rotor with a breathing crack model. The sectioned view of the cracked rotor, in both 
inertial and rotating coordinates, is shown in Fig. 1. The stiffness of the uncracked rotor system is symmetric (isotropic), 
and the damping due to the air resistance effect is assumed to be viscous. 


ft 1 



Figure 1. Section view of cracked Jeffcot rotor. 


The angular position of the unbalance vector can be expressed as 6(t) = cot + y/(l) + ft . where co is a constant spin 
speed of the shaft, y/(t) is a torsional angle, and fi is the fixed angle between the unbalance vector and the centerline 
of the transverse shaft surface crack. It should be noted that 

<S(t) = cot + y/(t), 0 = <i> = ® + (/>, and 0 = O =y/ (1) 


where <D is the spin angle of the rotor. 


The kinetic and potential energy for the rotor system subjected to lateral and torsional vibrations can be expressed as 
follows: 


T = -^JpO 2 +\M[i 2 + y 2 } + ^Ms 2 d 2 + Msd (■ -z sin 9 + y cos d) 

u = \{z y)*ify + \KV' 2 


where 


K/ = 


k k 

^ yz "'yy 


( 2 ) 

(3) 


(3a) 


is the shaft stiffness matrix in inertial coordinates, with the off-diagonal elements corresponding to the asymmetry 
introduced either by a crack. Using the Lagrange approach, the nonlinear, coupled equations for transverse and torsional 
motion for the rotor system take the following form: 


Mz + C/Z + k, z z + k zy y = F, +Ms[6 2 cos# + <9sin0j 
My + C,y + k v: z + k yv y = F v +Ms[d 2 sin 9 -9 cos 

J p y/ + C t if / + sC, (z sin 6-y cos 6)-s ( F. sin 0 - F y cos 6} +s ^[k- zz z + k zv y^ sin 6 -(k yz z + k vv y^j cos + 


dU 

dy/ 


(4) 


=m 


(5) 
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where F, and F y are the external forces (including gravity) in z and y directions, respectively, C, is torsional damping 
coefficient and T e is the externally applied torque. 

Uncracked shaft 

For the uncracked shaft the following conditions related to potential energy and stiffness components are satisfied: 


dU 

dy/ 


= k,y/ , 


k a =kyy=k and = k yz = 0 


Cracked shaft 

The stiffness matrix for a Jeffcott rotor with a cracked shaft in rotating coordinates can be written as: 



(k c 0) 


'k 


(A k c 0 

Kj? = 

0 k n 

\ *1 ) 

— 


,]-/(*) 

0 A k„ 

\ i J 


(6) 


( 7 ) 


where the first matrix refers to the stiffness of the uncracked shaft, and the second defines the variations in shaft 
stiffness A kc and A k in E, and /] directions, respectively. The function / ( <D ) is a crack steering function which 
depends on the angular position of the crack, <t> , and the selected crack model. 

The simplest crack model is the hinge model, where the crack is assumed to change from its closed to open state 
suddenly as the shaft rotates. The steering function for this model is defined as: 


/(®) 


JO for £<0 
[1 for £>0 


( 8 ) 


While the hinge model might be an appropriate representation for very small cracks, Mayes and Davies [8] proposed a 
model with a smooth transition between the opening and closing of the crack that is more adequate for larger cracks. In 
this case the crack steering function, or the Mayes modified function, takes the following form: 


/($) = 


l + cos(<J>) 


The stiffness matrix for a Jeffcott rotor with a cracked shaft in inertial coordinates, K y is derived as 


K, = TK^T = 


*zz Ky \ (k Oj /(<J>)k( / Ak 1 + Ak 2 cos2<t> Ak 2 sin2<t> 


yKz k yyJ 


0 k 


A k 2 sin 2<J> A/q - A k 2 cos 2<t> 


where the transformation matrix T is 


T = 


cos® -sin ® 


sin ® cos ® 


and 


A k x = - 


A k f + A k„ 


and A k 2 = - 


Ak c - A k„ 


( 9 ) 


( 10 ) 


(11) 


( 12 ) 


The cross-stiffness for deep cracks has been accounted for assuming Ak l? = Ak. j 6 . which based on relationship in Eq. 


7 A k 


5 


dU 


(12), yields A^j = — andAA: 2 = — . The derivative of the potential energy term in Eq. (5) takes the 


6 k 


6 k 


dy/ 


following form: 


5767-38 V. 3 (p.3 of 12) / Color: No / Format: Letter / Date: 2005-03-19 06:20:09 

3 

SPIE USE: DB Check, Prod Check, Notes: 


Please verify that (1 ) all pages are present, (2) all figures are acceptable, (3) all fonts and special characters are correct, and (4) all text and figures fit within the 
margin lines shown on this review document. Return to your MySPIE ToDo list and approve or disapprove this submission. 


dU__ _k_ qf(<$>) r + cos24>)z 1 + 2 zyAk 2 sin2<D + (A^ - cos2<t>) y 1 

ri// 4 8 i// L 

+ ./ { < 4>)k\k 1 |~^2 s j n 2<J) _ 2zy cos 2<t> - y 2 sin 2<1>J + k t yr 


(13) 


RESULTS AND DISCUSSION 

A steel simple supported shaft with a disc at the center is considered. A transverse crack is located near the disk. The 
first lateral critical speed is 240 rad/s. The parameters are shown in the Table 1. 

Table 1. Parameters for a simple-supported Jeffcott rotor with a crack. 


Disk mass, M 3 Kg 

Disk polar moment of inertia, J p 0.01 kg-m 2 

Eccentricity of the disk, s 2.2x10 ~ 5 m 

Shaft stiffness, k 1.728 x 10 5 N/m 

Damping ratio in lateral direction, 4) 

Damping ratio in torsional direction, 4< q 075^ 

Torsional natural frequency, co t 600 rad/s 


The equations of motion (Eqs. (4)-(5)) were numerically integrated using Runge-Kutta method with the time integration 
step of 2;r/(360o) . This time step was small enough to get an accurate solutions even for slow rotational speeds such 
as 25 rad/s (A t= 6.9><10" 4 s). Bifurcation diagrams and power spectra for both lateral and torsional responses were used 
for presentation and analysis of the nonlinear system dynamical behavior. The bifurcation diagram plots the rotor orbit’s 
z- (or y/-) coordinate (with a dot) for each shaft revolution as the keyphasor reference mark fixed on the rotor passes the 
same rotational angle. If the orbital motion were strictly synchronous, only the same “dot” would appear repeatedly. 



3 
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(a) (b) 

Figure 2. Bifurcation plots for cracked rotor: 4/ = 0.008 , = 0.0006 , Ak^/k = 0.4 ; (a) lateral response; (b) 

torsional response. 


Two bifurcation diagrams are presented in Figure 2, showing the calculated vertical and angular rotor displacements at 
each speed increment over the speed range from 5 to 260 rad/s, for the cracked rotor model with the crack’s depth of 
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0.4, and lateral and torsional damping ratios 4) = 0.008 and £, t = 0.0006 , respectively. The bifurcation plots appear to 
be alike; the results for the two different vibrational responses are comparable in terms of their general shapes. 
Resonances occur just below the lateral natural frequency and slightly above the half of lateral natural frequency, which 
one would expect for a non-linear system of this kind. Spanning the rotor spin speed in the range of 145-149 rad/sec and 
216-234 rad/sec, reveals the presence of a wide variety of the excited response characteristics. This is demonstrated in 
Figs. 3 and 4, showing the magnified “snapshots” of the sections of bifurcation diagram presented in Figure 2(a), and 
the corresponding orbital motions. During the chaotic behavior the lateral vibration amplitudes increase more than 1000 
times, which would be equivalent to the system failure. Figure 2(b) shows torsional response of the rotor which is not 
subjected to any externally applied torque excitation, and provides clear evidence for the existence of lateral/torsional 
coupling in the system due to presence of lateral excitations and crack. 



Figure 3. Bifurcation plot for the cracked rotor: = 0.008 , = 0.0006,Ak^/k = 0.4. 
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Figure 4. Bifurcation plot for the cracked rotor: = 0.008 , i; t = 0.0006 , Ak^/k = 0.4 . 


.00 


Figure 5 illustrates the lateral vibration responses for the uncracked and cracked rotor running at speed of 225 rad/s. 
There is no external torsional excitation applied to the rotor. It can be seen that uncracked rotor vibrates around its static 
deflection position with the small vibration amplitude (Fig. 5(a)). The corresponding power spectrum shows only 
synchronous frequency component due to unbalance excitation. In the presence of transverse crack of depth of 0.4, the 
rotor vibration amplitudes increase significantly, and the power spectrum illustrates the apparent system nonlinearity, 
indicating the multi-frequency content present in the vibration signal (Fig. 5(b)). 

Note the window of the FFT is taken to be the last steady state 50 system rotations. This condition was reached by 
allowing the inherent system hysteresis damp out all starting effects in the first 200 revolutions. If the motion is periodic 
or quasi-periodic, then the power spectrum will consist of a sequence of spikes at the fundamental frequencies, their 
harmonics, and the frequencies that are the sums and differences of the various frequencies. For the chaotic motion 
(e.g., non-periodic but non-random), the power spectrum shows a random broadband character. 

Figure 6 shows the torsional vibration responses for the uncracked and cracked rotor running at speed of 225 rad/s., 
without external torsional excitation applied to the rotor. Due to the nonlinear lateral/torsional coupling the unbalance 
excitation induces the uncracked rotor torsional response at predominantly torsional resonance frequency (600 rad/sec) 
and superharmonic frequency (2 at), (see Fig. 6(a)). After introducing transverse crack, the nonlinearities caused by the 
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breathing of the crack strengthen the modes coupling. There exist many frequencies in the torsional frequency domain. 
Unbalance excitation not only causes the amplitude of torsional vibration to increase 10 3 times larger than that without 
crack, but also significantly change the vibration patterns (Fig. 6 (b)). 




Figure 5. Lateral vibration response of rotor at speed 225 rad/s, = 0.008 , = 0.0006 ; (a) without crack; (b) 

with crack AkJk = 0.4 . 
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x to - 6 




Figure 6. Torsional vibration response of rotor at speed 225 rad/s, = 0.008 , = 0.0006 ; (a) without crack; (b) 


with crack Ak,/k = 0.4 . 


The comparison of lateral and torsional response for the cracked rotor is presented in Figure 7. Here, the values of 
parameters used for simulation are the same as those used for Figures 5 and 6, except for the increased value of 
damping, i.e., 4/ = 0.048 and 4) = 0.0754) • By comparing responses shown in Figs. 5(b) and 6(b) with the ones on 
Fig. 7, one can observe that the chaotic behavior of the response disappear for this values of the increased damping. 

Figures 8 and 9 illustrate bifurcation diagrams showing the response for different damping levels (4) = 0.005-0.05) 
and depth of the crack ( A kc Ik- 0.2 - 0.4 ), respectively. These results clearly suggest the possibility of the use of chaos 

mapping of monitored vibration signals as a diagnostic tool to detect some rotating machinery malfunctions. As the 
damping ratio or crack depth are ranged over the interval illustrated in Figs. 8 and 9, a wide variety of response 
characteristics are excited. 


Finally, a torsional excitation of T e (t) = 800sin(240?) was applied to the case of the cracked rotor, again, in addition to 
the unbalance and gravity forces. These results are shown in Fig. 10. The selected rotor spin speeds are an integer and 
not-integer fraction of the bending natural frequency, i.e., 0.1 co n =24 rad/s and 0. 1 3a> n =3 1 rad/s, respectively. For both 
speed cases the power spectra look similar, in addition to such crack “indicators” as the 2X and 3X components, the 
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Figure 7. Vibration of cracked rotor at speed 225 rad/s; = 0.048 , = 0.075^, , Ak»/k = 0.4 ; (a) lateral 


response; (b) torsional response. 
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Figure 8. Bifurcation plots for cracked rotor at speed to = 225 rad/s ; £, = 0.075^, AkJk = 0.4 ; (a) lateral 

response; (b) torsional response. 


5767-38 V. 3 (p.9 of 12) / Color: No / Format: Letter / Date: 2005-03-19 06:20:09 

9 

SPIE USE: DB Check, Prod Check, Notes: 


Please verify that (1 ) all pages are present, (2) all figures are acceptable, (3) all fonts and special characters are correct, and (4) all text and figures fit within the 
margin lines shown on this review document. Return to your MySPIE ToDo list and approve or disapprove this submission. 



(a) (b) 

Figure 9. Bifurcation plots for cracked rotor at speed io = 225 rad/s , = 0.01 ; % t = 0.075^,, Ak Jk = 0.4 ; (a) 

lateral response; (b) torsional response. 



xIO ' 5 



Frequency (rad/s) 




(a) (b) 

Figure 10. Orbits and power spectra of cracked rotor; ^ = 0.008 , \ t = 0.075^, Ak^/k = 0.4, 

r e (0 = 800sin(240t); (a) speed 24 rad/s; (b) speed 31 rad/s. 
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torsional excitation frequency co e = 240 rad/sec is seen along with the side frequencies co e ±co, co e ± 2co , and co e ±3co , 
which are located around co e . These frequencies are generated as a result of the non-linear mechanism of the breathing 
crack under the influence of unbalance, gravity, and harmonic torsional excitation. However, the analysis of the orbital 
motion reveals that in the case of speed of 31 rad/s the orbit rotates in the opposite direction of the rotor’s rotating 
direction. Simulations of orbits considering frequencies shown in corresponding to this case FFT spectrum (see Fig. 
10(b)) have shown presence of frequency components being backward whirl fraction. Therefore, for this speed, the rotor 
vibration in the presence of the applied external torque is not periodic but quasi-periodic, where the ratios of the 
involved frequencies are not the ratios of integers. 

CONCLUSIONS 

An analytical solution was developed for the Jeffcott rotor taking into account coupling between the torsional and 
bending vibrations. This was achieved by the inclusion of the third equation for torsional vibration. The modeling was 
further enhanced by considering the nonlinearities associated with a breathing crack located on the shaft adjacent to the 
disk. Numerical simulations of the uncracked and cracked cases were carried out by applying an unbalance excitation, 
gravity forces, and a torsional excitation. 

When comparing the torsional and lateral frequency spectrums for the uncracked and cracked cases, very clear 
distinctions were noticed. In each case (uncracked and cracked), the system was excited by the residual unbalance and 
gravitational forces, and then by the addition of a torsional excitation. The sum and difference frequencies have been 
observed in a response of cracked rotor due to unbalance and external torsional excitation. By employing nonlinear 
vibration theory, some phenomena evolving out of crack in rotating systems have been addressed, such as the numerical 
study of the bifurcation of the parameter dependent system. 

The uniqueness of the observed crack signatures, captured by the analytical approach, demonstrates an improvement in 
early diagnosis of rotor cracks. 

NOMENCLATURE 

C, Damping coefficients in z and y direction, respectively 

C, Torsional damping coefficient 

F , , F External forces in z and y directions, respectively 
J Disk polar moment of inertia 

M Disk mass 

co Constant spin speed of the shaft 

y/ Torsional angle 

P Angle between the unbalance vector and the crack center line 

k Stiffness of uncracked rotor. 

k , z , kyy Direct stiffness terms in z and y direction, respectively. 

k zy Cross stiffness term introduced by the shaft asymmetry 

k t Torsional stiffness 

s Eccentricity of the disk 

4/ Damping ratio in lateral direction 

4, Damping ratio in torsional direction 

co, Torsional natural frequency 

co n Lateral natural frequency 

T e (t ) External torque 

4, 4, H Rotor-fixed coordinates system 
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z, y Inertia coordinate system 

Me, Ak Reduced stiffness in £ and i) direction, respectively 
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ABSTRACT 

Critical components of propulsion systems frequently operate at high stress levels for long periods of time. The integrity 
of these parts must be proven by non-destructive evaluation (NDE) during various manufacturing steps and also during 
systematic overhaul inspections. Conventional NDE methods, however, have unacceptable limits. Some of these 
techniques are time-consuming and inconvenient for service aircraft testing. Impedance-based structural-health- 
monitoring (SHM) uses piezoelectric (PZT) patches that are bonded onto or embedded in a structure; each individual 
patch both actuates the surrounding structural area and senses the resulting structural response. The size of the excited 
area varies with the geometry and material composition of the structure. A series of experiments on simple geometry 
specimens (thin-gage aluminum square plates) was conducted for assessing the potential of E/M impedance method for 
structural damage detection. Based on the results of this preliminary study, further testing was conducted on a subscale 
disk specimen. Based on the results it can be concluded that the E/M impedance method has the potential to be used for 
damage detection of structures. The experimental method, signal processing, and damage detection algorithm should be 
tuned to the specific method used for structural interrogation. 

INTRODUCTION 

Structural health monitoring (SHM) is an essential element of modern structural components. The idea of equipping 
structures with sensors and actuators in an attempt to impart “smartness” has great potential for establishing a cost 
effective, in-situ maintenance routine 1 . This is especially true for critical high performance, aerospace components 
where accessibility using traditional nondestructive evaluation techniques is limited (e.g., turbine propulsion 
components) or even impossible (e.g., long duration space exploration vehicles). Some examples of manual and time 
consuming nondestructive evaluation techniques (NDE), which are typically utilized off-wing in the aerospace industry, 
include eddy current, ultrasonic, fluorescent penetrant, magnetic particle, optical, and radiographic inspection 2 . 
Regarding the cost effectiveness of reduced manual inspections, it is estimated that nearly 27% of an aircraft’s life cycle 
cost is spent on inspections and repairs 3 . With an on-line, self actuated system such costs can be dramatically reduced. 
Furthermore, the impact of such an in-situ SHM system is that it not only increases safety and performance, but also 
enables converting schedule based into condition based maintenance, thus reducing both down time and costs 4 . 

The impedance based approach has demonstrated unique features that meet the requirements of an on-line SHM system. 
This form of structural health monitoring is based on the use of a sensor/actuator patch to obtain real-time and 
continuous measurements that reflect the health status of the monitored structure. This paper describes the development 
of an impedance based SHM system from hardware acquisition, to software development, to procedural methodology 
and, lastly, experimentation using undamaged and damaged aluminum plates and disks. 



IMPEDANCE BASED STRUCTURAL HEALTH MONITORING 


Theory 

Impedance-based SHM uses piezoelectric (PZT : lead, zirconate, titarate) patches that are bonded onto or embedded in a 
structure. Each individual patch behaves as both an actuator of the surrounding structural area as well as a sensor of the 
structural response. The size of the excited area varies with the geometry and material composition of the structure. 
When a PZT material is subjected to an electric field it produces a mechanical strain, and when stressed it produces an 
electrical charge. For a PZT patch intimately bonded to a structure, driving the patch with a sinusoidal voltage sweep, 
for example, deforms and vibrates the structure. This is due to the patch applying a strain parallel to the structure’s 
surface. In reaction to these elastic wave inputs, the structure produces a localized dynamic response. This dynamic 
response is transferred back to the PZT patch, which sequentially produces an electrical response that is analyzed in 
regard to the impedance behavior 5 . The structure’s mechanical impedance is presented in the classical formulation as 


Z slr (co) = iojm e (co) + c e (co) - ik e (of)! co 


(1) 


The terms m e ,k e , and c, represent the structure’s mass, stiffness, and damping coefficients while co represents 
frequency. Due to mechanical coupling between the sensor and the host structure, this mechanical effect is picked up by 
the sensor and, through electro-mechanical coupling inside the active element, is reflected in the electrical impedance 
measured at the sensor’s terminals. Figure 1 illustrates the basic concept 6 . 
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Figure 1. Electro-mechanical coupling between the active sensor and the structure (Giurgiutiu et al. 1998). 


Solving the wave equation for the system shown in Figure 1 yields the total impedance as measured by the PZT sensor, 
Z (ajj , which contains both the structure’s Z slr ( (o) and the sensor’s Z P7T ( co) impedances: 


Z(co) = 


icoC 


1 - K 


31 


Zstr(M) 

Z s i r + Z pzT (CO) 


( 2 ) 


The term C denotes the zero-load capacitance, and /c 31 represents the coupling coefficient of the piezoelectric active 
sensor for in-plane vibration. At this point in time, the impedance based health monitoring technique is typically utilized 
as an empirical before-and-after tool that identifies changes in the damage state by noting phase shifts or magnitude 
alterations in the measured impedance as compared to a base line measurement. 



Experimental Set-up 

Figure 2 shows the experimental set-up utilized for capturing electromechanical (E/M) impedance data. The E/M 
impedance data were taken using an HP 4 1 94A impedance analyzer connected to a computer via the GPIB interface. A 
custom software program was developed for control of the impedance analyzer functions as well as acquisition, 
processing and display of the collected data. Of particular importance was the ability to modify the frequency sweep 
range, step size, excitation voltage, and equivalent circuit mode (series or parallel) for a given test. For each test, 
impedance data were acquired at every discrete excitation frequency within the sweep. In this study, only the real 
portion of the impedance was collected and analyzed as this has been shown to be the most sensitive in regard to 
damage identification 7 . Following data acquisition, the software provided functions for displaying collected data, 
averaging, performing damage metric calculations and exporting raw and reduced data. 


Computer 



Figure 2. Experimental setup. 


Damage Metric 

In order to quantify damage, a metric was employed that examines the change in the E/M impedance response with 
regard to the baseline condition. This damage metric, D, is expressed as 

D = {\~p ) 2 (3) 

where, 

_co v(XY)_ 

atX)a(Y) 


The term p represents the correlation coefficient concerning the baseline, (X). and damaged, (T), E/M impedance 
responses. It is a function of the covariance, cov(XY). and standard deviations, r>(X) and r>( Y). of the response values 8 . 
The damage metric will range from 0 for the undamaged case to 1 for fully damaged case (i.e., no correlation between 
the before-and-after impedance measurements). Similar approaches have been utilized in other studies (e.g., Giurgiutiu 
etal. 1998). 




EXPERIMENTAL PROCEDURE 


Plate Specimens 

A series of experiments on two, simple geometry specimens (thin-gage aluminum square plates) were conducted for 
assessing the potential of the impedance based health monitoring system. The plate specimens, each measuring 
lOOxlOO-mm with a thickness of 1.5-mm, were made from 6061-T6 aluminum alloy sheet. Each plate was instrumented 
with one lOxlO-mm PZT patch (material PSI-5A4E, thickness 0.19-mm attached using cyanoacrylate adhesive) 
purposely located at a vibration-sensitive location. This location was determined by conducting a finite element analysis 
focusing on the modal response of the plate and placing the patch in an area suspect to maximum displacements. After 
capturing baseline data for each of the undamaged plates, a single 10-mm straight, through-thickness EDM notch (0.3 
mm width) was utilized to simulate a crack. The notch locations for each plate are shown in Figure 3. As was the case 
for the patch placement, the notches were located in areas of maximum displacement as defined by the finite element 
analyses. During the experiments, the specimens were supported on foam to simulate free-free conditions. 

Preliminary tests were conducted to determine the frequency range best suited for measuring the damage dependent 
changes in the impedance response. Impedance values were collected for 15 repetitions over a wide, 20 - 200 kHz, 
frequency range for each specimen in the undamaged state. The 15 repetitions were then averaged (i.e., a single 
impedance plot using the 15 repeats) and compared to one of the 15 samples (randomly chosen) from within the same 
data set. This was done for each of the two plates. The intent was to look for a frequency window that produced the best 
repeatability over the 15 measurements in the undamaged state by avoiding noisy frequency regions that lack 
consistency between measurements. To this end, the damage metric was calculated for a 10 kHz wide, sliding window 
starting at 20 kHz and ending at 200 kHz. The step size for the window and the accompanying damage calculation was 
1 kHz. Once a suitable frequency range was selected, a second set of measurements utilizing this range were made on 
the plates following the introduction of the notch. 




Plate 1 Plate 2 

Figure 3. Schematic of specimens with simulated cracks and PZT sensors. 


Disk Specimen 

A second set of tests were conducted on a subscale disk that will be used to evaluate various propulsion system damage 
detection methodologies. The disk, shown in figure 4 was constructed of titanium with an outside diameter of 235 mm. 
and a bore diameter of 46 mm. The disk was supported in the same manner as the plate specimens. In this proof of 
concept evaluation, a 10mm square PZT patch was placed on the web of the disk. A set of 5 baseline impedance 
measurements were made over a range of 20-40 kHz. This range was selected using the same procedure as the plate 
samples. Following the baseline measurements, a small, 2.2 g., mass was attached the disk in order to simulate a 
change in the structural response of the disk. The mass was placed at 90 degree increments from the patch at the same 


distance from the center as the sensor and the impedance measurements repeated. Following the application of the 
mass, another set of measurements were made to determine the repeatability of the baseline E/M impedance response. 



Figure 4. Disk specimen with mass and PZT sensor locations. 


RESULTS 


Plate Specimens 

Figure 5 shows a plot of the damage metric versus the central frequency of the windowed region in the undamaged 
plates. As mentioned above, this procedure was done in order to search for frequency regions that offered good 
repeatability under constant conditions (i.e., undamaged case). It can be seen from this plot that the damage metric was 
consistently valued at zero from 20 kHz to approximately 45 kHz, meaning that the two plots correlate well in that 
frequency range. The damage metric, in regions of good repeatability, should be close to zero, since the material state 
was maintained (zero damage in this case). Beyond 45 kHz, however, there was a measurable increase in the damage 
parameter. Based on these results, a frequency range of 20-40 kHz was chosen for subsequent E/M impedance 
measurements. 



Figure 5. Windowed correlation coefficient (20-200 kHz, 10 kHz window, 1 kHz step) for two uncracked plates. 

Next, a comparison was made between the baseline of plate 1 and the baseline of plate 2. Figure 6 shows the E/M 
impedance plots as well as indicating the damage parameter value based on the correlation of the two plates (Equation 
3). Each of the impedance plots was an average of 15 repetitions. The damage metric (.0=0.965) indicated that there was 
much discrepancy between the two baseline measurements. As stated earlier, the large specimen to specimen variations 
dictate that the method be utilized as a comparison of baseline and subsequent measurements within the same specimen. 
Similar observations were made by Winston et al. 200 1 2 . 



Figure 6. Comparison of undamaged plates 1 and 2 (20-40 kHz). Damage metric, D, was 0.965 indicating a lack 

of correlation between the two plates. 




Figure 7 shows the results of comparing the baseline and the notched (i.e., damaged) cases for plate 1. The presence of 
the notch did modify the frequency response function. This was indicated by the shifting of the resonant frequencies and 
by the appearance of new resonances. The calculated correlation coefficient for this case was jO=0.497, and the damage 
metric was D= 0.253 based on Equation (3). Next, Figure 8 shows the results for plate 2 comparing the undamaged 
baseline with the damaged case, again, calculated using 15 repetitions for each plot. The correlation coefficient for plate 
2 was /i=0.365, and the damage metric was D= 0.403. Note that plate 2 before-and-after notch comparisons showed a 
larger damage metric as compared to the plate 1 case, even though the plate 2 notch was located farther from the PZT 
patch. At this point there is no explanation for this behavior, although, the uniqueness of the plate to plate behaviors 
probably plays a part. Furthermore, the changes in the E/M impedance could also be influenced by the notch location as 
it relates to the plate’s vibration displacement response and node locations. Further research needs to be conducted 
concerning these results. 



Figure 7. Averaged E/M impedance results for plate 1 with and without a notch. 



Frequency (kHz) 


Figure 8. Averaged E/M impedance results for plate 2 with and without a notch. 




Disk Specimen 

Figure 9 shows a comparison between the baseline E/M impedance response of the disk and the disk with the mass 
applied at one of previously specified locations. Examination of the response plot shows both changes in amplitude and 
shifting of response peaks. Figure 10 displays a summary of the calculated damage parameters at the various locations 
compared to the baseline condition. Results show that application of the added mass produces a measurable change in 
the damage parameter at each location. Measurement of the E/M response after testing indicates that the system returns 
to the baseline values (i.e. a zero damage parameter). Based on these results, it seems that the damage parameter is 
smaller (indicates less damage) with the mass at the farthest location. This appears to indicate that the location and 
number of sensors in an application must be considered in choosing a health monitoring scheme in order to assure full 
coverage of the structure with adequate sensitivity. 
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Figure 9. Averaged E/M impedance results for disk with and without added mass. 



Figure 10. Damage metric with added mass at various locations. 

CONCLUSIONS 

The results of this study demonstrate the ability of impedance based health monitoring to detect damage and structural 
changes in both aluminum plates and subscale propulsion type components. In addition, a methodology was developed 
that assists in choosing a frequency range for E/M impedance measurements in order to reduce test time and increase 




data integrity. This method was based on finding frequency ranges that offered the best repeatability (i.e., lowest noise) 
as expressed by either a high test to test correlation or low damage value. Next, simulated cracks, in the form of narrow 
notches, were analyzed using PZT patches as impedance measuring transducers. Conclusions from the experiments on 
the two aluminum plates were as follows: The assessment of the impedance behavior regarding the two undamaged 
plates indicated that variations between specimens overrode changes due to damage. Therefore, it was shown that the 
impedance method was best utilized as a comparison of baseline and subsequent measurements concerning the same 
specimen. Next, the E/M impedance measurements were able to recognize the changes between the undamaged and 
damaged (i.e., notched) cases for each plate, although, larger changes were observed for plate 2 even though the notch 
was located farther from the PZT patch as compared to plate 1 . This may have been due the uniqueness of the plate to 
plate behaviors, as well as the influence of notch and patch locations as it relates to the plate’s displacement response 
and node locations. Further research is needed regarding this observation. Based on the results of the plate tests, the 
method was applied to a subscale turbine disk. The results from these tests indicate that the method described above has 
potential application in the structural health monitoring of propulsion system components. Although the tests were 
conducted on an idealized component in a static condition, it is reasonable to assume that, based on the demonstrated 
sensitivity, the method can be extended to larger scale components subjected to dynamic loads. Due to the sensitivity of 
the measurements, however, it may be necessary to consider utilizing an unbounded damage parameter such as the one 
employed by Peairs 5 . 

With the establishment of the hardware and software as well as success concerning the preliminary experiments, the 
system will be further tuned to the various material systems (e.g., composites or thermal protection system materials) 
and structures that are of most interest to NASA. This includes the application of the technique to rotor components 
(either aero or space propulsion) by utilizing slip ring connections or interrupted testing. In addition, in cooperation with 
the Virginia Polytechnic Institute and State University, research is being conducted regarding wireless and self 
powering applications especially relating to space exploration vehicles. This will allow impedance based SHM to be 
applied in a vast array of applications. Lastly, the same PZT patches can be utilized as ultrasonic transducers, thereby 
creating a layered health monitoring system that is permanently attached to a structure allowing for continuous feedback 
regarding system health. This will allow for the identification of a damage condition and, based on ultrasonic response, 
identify the size, location and type of damage. 
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ABSTRACT 

The coupling between lateral and torsional vibrations has 
been investigated for a rotor dynamic system with breathing 
crack model. The stiffness matrix has been developed for the 
shaft element which accounts for the effect of the crack and all 
six degrees of freedom per node. Since the off-diagonal terms 
of the stiffness matrix represent the coupling of the respective 
modes, the special attention has been paid on accurate 
determination of their values. Based on the concepts of fracture 
mechanics, the variation of the stiffness matrix over the full 
shaft revolution is represented by the truncated cosine series 
where the fitting coefficient matrices are extracted from the 
stiffness matrices of the cracked shaft for a number of its 
different angular positions. The variation of the system 
eigenfrequencies and dynamic response of the rotor with two 
cracks have been studied for various shaft geometries, crack 
axial locations, and relative phase of cracks. 

INTRODUCTION 

Any new advancement in on-line detection and diagnosis 
of critical malfunctions in rotating machinery can be extremely 
beneficial to industry [1]. This especially applies to fatigue 
cracks on the shaft, which present a potential source for 
catastrophic failures of rotating machinery. Most reported 
studies are focused on two crack signatures, i.e., twice the 
running frequency component (2X) and the subharmonic 
component at approximately half of the shaft critical speed. A 
study done by Tondl [2] was one of the first to investigate the 
effect of the coupled lateral and torsional vibrations on 
turbogenerator rotor stability. He concluded that the combined 
effect of torsional stiffness and bending stiffness results in 
speed intervals where the rotor vibrations due to residual 
unbalance become unstable. A comprehensive literature review 


Jerzy T. Sawicki 

Rotor-Bearing Dynamics & Diagnostics Laboratory 
Cleveland State University 
Cleveland, OH 44115-2214 

George Y. Baaklini 

NASA Glenn Research Center 
21000 Brookpark Road, MS 77-1 
Cleveland, OH 44135 


of various crack modeling techniques and system behavior of 
cracked rotor is given by Wauer [3]. Papadopoulos and 
Dimarogonas [4] used a non-rotating cracked Timoshenko shaft 
to demonstrate the existence of an apparent coupling of 
torsional and bending vibration. They modeled crack using the 
local flexibility matrix, and then proceeded to study the 
vibration spectra in the presence of harmonic excitations. Other 
researchers also dealt with the problem of coupled vibrations, 
for example Plaut and Wauer [5] investigated resonances and 
instabilities in coupled flexural and torsional vibrations of a 
rotating shaft. Muszynska et al. [6] analytically and 
experimentally analyzed lateral/torsional coupling mechanisms 
resulting from combinations of unbalance, shaft stiffness 
asymmetry, and radial sideload. Muszynska experimentally 
observed torsional resonance at speeds equal to 1/8, 1/6, 1/4, 
and 1/2 of the lowest torsional natural frequency. Bently et al. 
[7] continued this study with special attention paid to the 
analysis of a “snapping” action which occurs when during 
rotation the rotor experiences a peak torsional acceleration. 

The topic of cracked rotor vibrations has been analyzed in 
a number of published works [8-16]. They have been focused 
on the study of dynamic behavior of rotors with the so-called 
breathing type of crack during the passage through the critical 
speed at the constant angular acceleration or deceleration. For 
example, Sawicki et al. [9, 10] studied the accelerating cracked 
rotor response using the angle between the crack centerline and 
the rotor whirl vector to determine the closing and opening of 
the crack. This allows one to study the rotor dynamic response 
with or without the rotor weight dominance by taking into 
account nonsynchronous whirl. Henry [11] investigated how 
the gravity and the unbalance affect the vibration response of a 
cracked shaft. A theoretical and experimental study of the 
effects of a transverse crack on the rotor dynamic system was 
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given by Mayes [12]. From the experimental results, it is 
observed that the crack has most significant effects on the 
response of the rotor dynamic system when the phase angle 
between imbalance and crack is in the range of 45° and 135°. 
Out of this range, the rotor behaves like an un-cracked rotor. 
Gash [13, 14] provided a comprehensive investigation of the 
stability behavior and insightful study of harmonic resonances 
of a cracked Jeffcott rotor with the hinge model of the crack. 
Mayes and Davies [15] modified hinge model to account for 
deep cracks, by introducing the crack cross flexibility. Collins 
et al. [16] used axial impulses for crack detection in rotating 
shafts. 

Practical simulation of cracked rotordynamic systems calls 
for application of finite element analysis to account for 
geometric characteristics of the rotor. Papadopoulos and 
Dimarogonas [17] derived flexibility matrix for the shaft 
element with open crack. Later, they studied [18] coupling 
between bending, longitudinal and torsional modes of vibration 
for non-rotating shaft with an open crack. Sekhar [19] 
presented results of finite element analysis of the flexural 
vibration response of the cracked rotor with two open 
transverse cracks, focusing mainly on the stability study and 
eigenfrequency analysis. Flowever, an open crack model in a 
case of rotating shaft is not practical. The rotor vibration 
characteristics with such crack model can be very different than 
with the breathing crack model. Darpe et al. [20] presented 
study of coupled longitudinal, lateral, and torsional vibrations 
for the cracked rotating shaft using a response-dependent non- 
linear breathing crack model. The signs of the overall stress 
intensity factors (SIF) at any point along the crack edge are 
used to judge whether the crack is open or close, thus 
determining the crack closure line position. 

In this paper, the coupled lateral and torsional vibrations of 
a rotating shaft with two breathing cracks are investigated using 
the finite element approach. Based on the concepts of fracture 
mechanics, the variation of the stiffness matrix over one rotor 
revolution is expressed by the truncated cosine series, in which 
the fitting coefficient matrices are determined from the stiffness 
matrices of the cracked shaft at five different angular positions. 
The method accurately predicts not only the direct stiffness 
terms, but also the off-diagonal terms which account for 
coupling mechanisms of the respective modes. The developed 
approach accurately describes crack breathing action and can 
be effectively applied to the frequency analysis of the rotor 
with two cracks at two different angular and/or axial positions, 
and for different shaft geometric ratios. Coupled lateral and 
torsional vibrations of a rotating shaft with two breathing 
cracks of different relative phase, under unbalance, gravity and 
external torque excitations are studied. 

NOMENCLATURE 

a Depth of the crack 

a Normalized depth of the crack ( a / R ) 
b Flalf-width of the crack 

b Normalized half-width of the crack ^2 a - a 2 

c c Cracked element flexibility matrix 
C Global damping matrix 

D Shaft diameter 

E Modulus of elasticity 

F Global force vector 


G Shear modulus of elasticity 
I 6x6 Identity matrix 

J Strain Energy Density Function 
J P Polar moment of inertia 
J D Diametral moment of inertia 
k c , K ( Cracked element stiffness matrix in rotor-fixed 
and inertial coordinate systems, respectively 
K ( cot ) Global stiffness matrix 

/ Length of shaft element 

L Length of rotor 

M Disk mass 

M Global mass matrix 

q Displacement nodal vector 

R Shaft radius 

t Instantaneous time 

T e Torsional excitation 

k Shape coefficient for circular cross section 

v Poisson’s ratio 

p Material density 

s Unbalance eccentricity 

co Shaft spinning speed 

co e External excitation frequency 

6 Rotational angle of the rotor 

x, rj, c Rotor- fixed rotating coordinate system 

x, y, z Inertial coordinate system 

MODEL OF THE CRACKED SHAFT ELEMENT 

Figure 1 shows a shaft element with a transverse crack of 
depth a, at distance x from node 1, loaded with axial forces 
P\, P-i , shear forces P 2 . P$ , P H ■ P t , bending moments F\. P 6 , 
P\ i > Pn> and torsional moments P 4 ’PlO ■ All six degrees of 
freedom per node are considered. 



Figure 1. Cracked shaft element in general loading. 

Since the strain energy induced at the tip of the loaded 
crack introduces the considerable local flexibility in the specific 
cracked beam element, the crack affects only the stiffness 
matrix of the given finite element. The geometry of shaft 
section with partially open crack is shown in Fig. 2, where a is 
the depth of the fully open crack and b is the half-width of the 
crack. 

The node displacement in the direction of load P t , induced 
only due to the presence of the crack with depth a, can be 
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calculated as follows [19]: 



Figure 2. Geometry of the shaft section with partially open 
crack. 


where A is the area of the cracked section of the shaft (see 
Figure 2) and J is the Strain Energy Density Function expressed 
as [21]: 


J(A) = 


(6 } 

I* 

V i = 1 


T 


( 6 V 


V i=i y 


+ (l + v) Y,K m 


V i=l 


(2) 


Here, E is the modulus of elasticity, v is the Poisson ratio, Ky 

are the crack Stress Intensity Factors (SIF) for i = I, II, III crack 
displacement modes and load index j = 1,2,. . ., 6 . 

For a strip of unit thickness with a transverse crack the 
values of SIF are known [19] and thus, employing integration 
along the tip of the crack with the variable crack depth, yields 
the approximate values for SIF in Eq. (2). Based on Eq. (1), the 
local flexibility due to the crack is calculated as [22]: 

(3) 

The expressions for terms of the total flexibility matrix of the 
cracked element c c = [ c i/] (ij =1,2,.. .,6) are given in 
Appendix A. 

MODEL OF THE CRACK BREATHING ACTION 

For the majority of rotating machinery, the static deflection 
is much larger than the rotor’s vibration amplitude. Under this 
assumption, the instantaneous shaft angular position can be 
employed to determine the amount of crack opening. Two 
extreme scenarios exist when the crack is either closed (0 = 0) 
or open ( 9 = n ). For the spinning rotor at constant speed, a 
truly breathing crack behavior is represented by gradual 
opening of the crack from its fully closed to the fully open state 
and visa-versa. Since the flexibility of the cracked shaft section 
changes with different amount of crack opening, the concept of 
Crack Transition Line (CTL) is introduced. The CTL, which 
separates the open and closed portions of the crack, as 
illustrated in Fig. 2, is an imaginary line perpendicular to the 
crack edge. The integration limits for evaluation of flexibility 
coefficients depend on the amount of crack opening. 


The crack edge is divided into (ATI) points. For any shaft 
angular position 0 = i |i80°/a| , the CTL is located at 

b-2ib/N (i = 0...A) along the crack edge. The terms of 
flexibility matrix are continuously updated for each angle 
<9 = 180 °/a . In this way CTL keeps changing with different 
angles of shaft rotation, and thus reflecting the true crack 
“breathing” action. 

The stiffness matrix for the shaft cracked element in the 
rotor-fixed coordinates system is calculated as: 

[ k mn\ =[ c ™], * m,n = 1>2,...,6; i = 0...A (4) 

where formulas for [c mn ] are given in Appendix A. 

In order to ensure the effectiveness of simulation, the 
variation of the cracked element stiffness with time (angular 
position), can be expressed by the truncated cosine series [23]: 

4 

k c = k 0 t ^k ( cos(iot) (5) 

1=1 

The coefficient stiffness matrices ky, j = 0,1,. ..4 , are 

calculated based on the assumed amount of the crack’s opening 
at its specific (five) angular positions. 

Let define that stiffness matrices k i(c , k , k 2 , k 3l? ,and 

k op to correspond to the rotor angular location at 0, n/4 , 7t/ 2 , 
3^/4 , and n , respectively. They are calculated by assigning 
value of i' = 0, A/4, A/2, 3A/4 , and A, respectively, in 
Eqs. (A2), i.e., by solving the following set of equations: 

k c (Jlt = 0 — k HC 5 k c I K — k tf ’ k c | n — k 2 
COt=— (Ot=— 

4 2 

k c 1 3»- = k 3? , k c \ wl=rr = k (6) 

cot = — 

4 

As a result, the fitting coefficient matrices take the following 
form: 

. k 2? k k k k 3 

k n = — — h — — t — - + — + — - 
4 8 8 4 4 

k i = ^( k » c - k «p + >/2k 9 - V2k 3? ) 
k 2 =^(- 2k 2 9 + k „ f + k op) (7) 

k 3 = ^( k «c - k op - ^ k , + V2k 3? ) 

k _ k 2 ? | k i)C | k op k ? k 3? 

4 4 8 8 4 4 

The developed method predicts very well the variations of 
all terms of cracked element stiffness matrix due to the 
“breathing” action of the crack, over one shaft revolution (see 
Fig. 3 1 ). Papadopoulos and Dimarogonas [17, 23] proposed 
similar approach but they assumed that the stiffness terms reach 
minimum values for the fully open crack, which is correct for 


1 Figure 3 is shown in Appendix B. 


3 


Copyright © 2005 by ASME 


diagonal coefficients ( k n , i = 1,...,6) but not for all cross- 
coupled terms. Some of the cross-coupled stiffness coefficients 
reach maximum or minimum values when the crack is half- 
open or half-closed, and become zero when the crack is fully 
open (see Fig. 3). 

From the finite element static equilibrium conditions, 12 
degrees of freedom of one element can be written as 

(q l ,q 2 ,-,qn) T = ( 8 ) 

where 
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The shaft is divided into 8 equal-length finite elements and two 
shown cracks shown are located at element 6 and 7, near the 
right disk (see Fig. 4). 


Table 1: Numerical Model Physical Parameters. 



Physical parameter 

Value 

Units 

L 

Shaft length 

1.12 

m 

D 

Shaft diameter 

0.03 

m 

p 

Material density 

7750 

Kg/m 3 

V 

Poisson’s ratio 

0.3 


E 

Modulus of elasticity 

2.07x10" 

N/m 2 

G 

Shear modulus of elasticity 

7.96xl0 10 

N/m 2 

M 

Disk mass 

3 

kg 

Jp 

Disk polar moment of inertia 

0.018 

kg m 2 

Jd 

Diametral moment of inertia 

0.01 

kg m 2 

s 

Unbalance eccentricity 

5.4xl0" 5 

m 


Next, the stiffness of the cracked element in inertial coordinates 
system can be found as 


K„ = TTjk TT t 


(10) 


where T e is the transformation matrix from rotor-fixed 
coordinates system to the inertial coordinates system 


and: 


T„ = 


f T 

a 

o 


V 



"l 0 0 N 

0 cos(ot) sin(®t) 

0 -sin(fttf) cos(ot) y 


( 11 ) 


( 12 ) 


When assembling the global stiffness matrix for the shaft, 
the cracked element stiffness matrix K r replaces the stiffness 
matrix of the corresponding uncracked element. The equations 
of motion for the complete rotor system in an inertial 
coordinates system take the following form: 

Mq + Cq + K(ot)q = F (13) 


where M, C, and K {cot) are the mass, damping and stiffness 
matrices, respectively, for the whole rotor system. The stiffness 
matrix is continuously updated with the angular position cot of 
the shaft. The force vector F can contain any kind of forces 
and moments applied at any nodes in the global inertial 
coordinate system. 

NUMERICAL SIMULATIONS 

A two-disk rotor system considered for numerical studies 
is shown in Fig. 4. Physical parameters of the model are listed 
in Table 1. The boundary conditions for the model follow the 
same conditions as for simply supported beam except constrain 
imposed on torsional degree of freedom at the right support. 


crack 2 


I 


LL 

/ 


A 


0.896 m 


1.12m 


Figure 4. A two-disk rotor system with two cracks near the 
right disk. 

Frequency Analysis 

Frequency analysis is carried out for an undamped non- 
rotating rotor. Three cases are considered for the cracks located 
at different axial locations, different phases, depth, and shaft 
geometric ratios. 

Case 1 : 

Two open cracks of depths of a l /D = a 2 /D = 0. 4; the 
position of one crack is fixed (0.756 m from the left support, 
see Fig. 4) while the position of the other one keeps changing 
from the left to right end. 

The variation of the first two mode frequencies, as a 
function of axial position of the second crack, is illustrated in 
Fig. 5(a). When the second crack is near the ends of the rotor 
the frequencies are almost the same as for the rotor with only 
one crack. The effect of the second crack on the rotor first two 
mode frequencies gradually increases, reaching the maximum 
as the crack’s position approaches the middle of the shaft. 

The third mode frequency, shown in Fig. 5(b), 
predominantly torsional, is very little affected by the position of 
the second crack. One can observe that when this position is 
between the left support and the left disk, the mode frequency 
is almost identical as for the rotor with one crack, i.e., 598 
rad/s. For the position anywhere between two disks, the 
frequency drops to 595 rad/s, and for the position between right 
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disk and the right support the frequency is reduced further to 
591.5 rad/s. 



(a) 


600 


598 

■3 

ro 

o 596 

Q) 

cr 

0) 

^ 594 - 

■D 

592 

0 0.2 0.4 0.6 0.8 1 1.12 

Second crack position (m) 

(b) 

Figure 5. Variation of mode frequencies for various 
positions of the second crack: (a) first two mode 
frequencies, and (b) third mode frequency. 




Second crack phase (deg.) 

(b) 

Figure 6. Variation of mode frequencies for various values 
of phase and depth of second crack: (a) first mode 
frequency, and (b) third mode frequency. 


Case 2: 

Two cracks have fixed axial positions; one crack at the 
distance 0.756 m and the second one at 0.896 m, measured 
from the left support (see Fig. 4). The first crack, of fixed depth 
a l /D = 0.4 , is always open, while the second is allowed to be 
oriented at any angle (phase) with respect to the positive z-axis, 
and have different depths. The changes of rotor mode 
frequencies with different depths and angular position of the 
second crack are illustrated in Fig. 6. It should be noted that the 
crack center line is defined as a line perpendicular to the crack 
edge, and here, the crack phase angle is defined as the angle 
between the positive z-axis and the crack’s centerline (see Fig. 
2). As expected, that the maximum reduction of frequencies 
happens to be when both cracks are fully open (phase 180°). 
Also, with the growing depth of the second crack, the rotor 
frequency reduction becomes more significant. 


Case 3 : 

Two cracks are located at fixed axial positions. Both cracks 
have fixed depths ( a 1 /D = a 2 /D = 0.4 ), but the shaft has 
different geometric ratios, i.e., for the given length of the shaft 
( L = 1.12 m) the shaft diameter is varied ( D = 0.2, 0.08, 0.06, 
0.04, 0.02 m). Also, the phase of the second crack keeps 
changing. As can be seen on Fig. 7(a), the third, predominantly 
torsional mode frequency is not much affected by the crack’s 
phase. To emphasize the scale of possible changes Figs. 7(b) 
and 7(c) show changes of the normalized frequency (with 
respect to its maximum value). For the rotor with smaller 
geometric ratios (i.e., larger shaft diameter), the variation of the 
frequencies is more sensitive to the phase change of the crack. 
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Figure 7. Variation of mode frequencies for different shaft 
geometric ratios: (a) variation of the third mode frequency 

(b) variation of the normalized frequency for the third mode, 

(c) variation of the normalized frequency for the first mode. 


Lateral and Torsional Vibration Response 

In this part of investigation two breathing cracks (rotating 
rotor) located at element 6 and 7 have the same depth 


Oj / D = a 1 / D = 0.4 . The stiffness variation property of the 
cracked element has been described in the previous section. 
The bending and torsional frequencies of the uncracked rotor 
are 214.4 rad/s and 601.3 rad/s, respectively. 

In order to explore the crack-induced mode coupling 
phenomenon, an external sinusoidal excitation torque 
T e = 600sin(<y e t) is applied at the first disk (node 3), with 
w e = 322 rad/s as an excitation frequency. The selected shaft 
spinning speed is a> = 43 rad/s. The purpose is to investigate 
how the torsional excitation affects the lateral vibration 
response of cracked rotor. Vibration responses for the three 
different phases of the two cracks are shown in Fig. 8 2 . 

In all these plots of Fig. 8, the existence of the first 
harmonic (IX) and super-harmonics (2X, 3X, 4X,) of the shaft 
running frequency in lateral frequency spectrum provide strong 
indication of the cracks’ presence. One can also notice the 
existence of torsional frequency, co e . in the lateral frequency 
spectrum, which clearly demonstrates the coupling mechanism 
between the lateral and torsional vibration modes. The 
appearance of sum and difference frequencies ( co e ± na> ) 
around the torsional excitation frequency is the result of 
interaction of the torsional excitation and the synchronous 
frequency and its higher harmonics (2X, 3X, 4X). 

For the case with cracks’ phase difference of 0°, shown in 
Fig. 8(a) and 8(b), the two cracks simultaneously open or close 
during the rotor’s rotation.. In this case the vibration spectrum 
of the two-crack rotor is similar to the one-crack rotor’s 
behavior, except the increased amplitude due to higher stiffness 
reduction. 

For the case with phase difference of 90°, shown in Fig. 
8(c) and 8(d), whenever one crack is closed, the other one is 
always half-opened and half-closed. While 2X, 3X and 5X 
peaks remain about the same, as in Fig. 8(a) and 8(b), the 
amplitude of 4X peak is dramatically decreased in the 
frequency spectra, and the vibration amplitude is slightly 
reduced. 

The third case with the phase difference of 180° implies 
that whenever one crack is fully opened, the other one must be 
closed. In this case the vibration amplitudes are reduced further 
than in the previous case. Although the super-harmonic 2X and 
4X peaks are similar as in Fig. 8(a) and 8(b), IX, 3X and 5X 
peaks are appreciable decreased to their minimum values. 

Due to the lateral and torsional coupling mechanism 
induced by the presence of cracks, the externally applied 
torsional excitation significantly affects the rotor lateral 
vibration behavior. This is clearly shown in Figs. 9 3 (a) and (b), 
where the vibration response orbits are presented without and 
with torque, respectively. The Fig. 9(b) shows the orbit that 
rotates (thin line) in the opposite direction of the rotor’s 
rotating direction. Simulations of orbits considering frequencies 
shown in corresponding to this case FFT spectrum (see Fig. 
8(b)) have shown presence of frequency components being 
backward whirl fraction. The vibration in the presence of the 
applied external torque is not periodic but quasi-periodic, where 
the ratios of the involved frequencies are not the ratios of 
integers. 


2 Figure 8 is shown in Appendix B. 

3 Figure 9 is shown in Appendix B. 
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CONCLUSIONS 

Coupled lateral-torsional vibration finite element analysis 
of a rotor with two breathing cracks has been conducted for 
frequency analysis and dynamic response to excitation forces 
such as mass unbalance, weight, and the external torque. The 
multi-crack rotors have much more complicated stiffness 
asymmetry behavior being parametric, time-varying systems 
with dissimilar inertia moments in two perpendicular directions 
and at various axial shaft locations. The issue of accurate 
modeling of the stiffness matrix with the periodic coefficients 
in the global inertial coordinates system is crucial to account 
for the proper breathing behavior of multi cracks located at 
various axial and angular positions. 

The developed model for the “breathing” cracked element 
is based on fracture mechanics. The direct stiffnesses, as well 
as cross coupling stiffnesses are estimated as the crack opens 
and closes. The stiffness matrix variation due to the crack 
breathing is determined by the stiffness matrices of cracked 
shaft with the crack located at five different angular positions. 
Such an approach simplifies frequency and dynamic response 
analysis regardless of a number of finite elements and cracks’ 
location. 

The axial location of cracks affects the rotor frequencies. 
Once the two cracks are located close to the middle section of 
the shaft, the maximum reduction of frequencies occurs. The 
frequencies also depend on the relative orientation of two 
cracks and the shaft geometric ratios. The variation of rotor 
mode frequencies is the most significant for the small shaft 
geometric ratios. The modes which affected the most are 
predominantly bending modes. If the two breathing cracks are 
in phase, i.e., both of them are simultaneously either opened or 
closed, the vibration behavior pattern is almost the same as for 
the rotor with one crack, except the change in vibration 
amplitude. For the case of two out-of-phase cracks, the 
vibration spectrum presents different signatures, with 
diminished IX, 3X and 5X harmonics. Finally, for the cracks’ 
phase difference of 90°, the peak of super-harmonic frequency 
4X is dramatically reduced in the vibration spectrum. It is 
noteworthy that the 2X peak stays the same all the time. 

The orbit analysis reveals that there are specific crack 
signatures due to the vibration mode coupling. When the 
external torque is applied the mode coupling induced by the 
crack causes vibration to become quasi-periodic, with the 
response orbit rotating in the direction opposite to the direction 
of rotation of the rotor. 

The presented method is applicable to analyze the dynamic 
response of rotating shafts with two or more breathing cracks. 
The uniqueness of the observed crack signatures, captured by 
the analytical approach, presents potential in early diagnosis of 
rotor cracks. 
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Figure 3. Variation of stiffness coefficients over one revolution of a cracked shaft element 
(a/D = 0.4, 1 = 0.14m, D=0.03m). 



(a). Phase difference 0° ; no torsional excitation. (b). Phase difference 0° ; torsional excitation 
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(c). Phase difference 90°; no torsional excitation. 



(d). Phase difference 90°; torsional excitation. 




(e). Phase difference 180°; no torsional excitation. (f). Phase difference 180°; torsional excitation. 

Figure 8. Vibration response FFT of a rotor with two breathing cracks of different phase. 
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Figure 9. Vibration response orbits with phase difference 0°: (a) without torsional excitation, (b) with torsional excitation. 
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